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ABSTRACT 

We have acquired sub-millimeter observations of 33 fields containing 37 Herbig 
Ae/Be (HAEBE) stars or potential HAEBE stars, including SCUBA maps of all 
but two of these stars. Nine target stars show extended dust emission. The 
other 18 are unresolved, suggesting that the dust envelopes or disks around these 
stars are less than a few arcseconds in angular size. In several cases we find that 
the strongest sub-millimeter emission originates from younger, heavily embedded 
sources rather than from the HAEBE star, which means that previous models 
must be viewed with caution. These new data, in combination with far-infrared 
flux measurements available in the literature, yield SEDs from far-infrared to 
millimeter wavelengths for all the observed objects. Isothermal fits to these 
SEDs demonstrate excellent fits, in most cases, to the flux densities longward 
of 100 /im. We find that a smaller proportion of B-type stars than A and F- 
type stars are surrounded by circumstellar disks, suggesting that disks around 
B stars dissipate on shorter time scales than those around later spectral types. 
Our models also reveal that the mass of the circumstellar material and the value 
of f3 are correlated, with low masses corresponding to low values of 0. Since low 
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values of (3 imply large grain sizes, our results suggest that a large fraction of the 
mass in low-beta sources is locked up in very large grains. Several of the isolated 
HAEBE stars have disks with very flat sub-millimeter SEDs. These disks may 
be on the verge of forming planetary systems. 

Subject headings: ISM: clouds - (Stars:) circumstellar matter - Stars:formation 
- Stars: pre- main sequence - Stars: variables: T Tauri, Herbig Ae/Be - Submil- 
limeter: stars 



1. Introduction 



sequence stars ( 


Herbidll960t 


Strom et al. 


1972: 


1992; 


The, de Winter and Perez 


1994 


). Since 



)er- 



ship in this class, our understanding and consequently our classification rules for HAEBE 
stars have changed somewhat to include some stars of slightly later spectral class and oth- 
ers that do not illuminate reflection nebulae, including "isolated" HAEBE stars. At visible 
wavelengths, spectra of HAEBE stars often include broad emission lines that show rapid vari- 
ability while at longer wavelengths HAEBE stars often are associated with infrared and far- 
infrared excesses and millimeter wavelength line emission. These properties are very strongly 
associated with circumstellar gas and dust at a wide range of temperatures. Most workers in 
the field now accept that the bulk of the circumstellar gas and dust aroun d many, if not most 
of th e se stars, is found in fairly massive (^O . Ol Mm) circumste ll ar di sks (H ami douche et al. 



2006; Eisner et al. 



2003 



Meeus et all l200lt iNatta et all l200lt iNatta. Grinin. fc Mannings 



2000l ; Hillenbrand et al.lll992l ). HAEBE stars are the progenitors of Vega like stars. Prob- 
ing the circumstellar disks of HAEBE stars is important to investigate the likelihood of 
planet formation around intermediate mass stars. Accretion disks around the more massive 
Herbig Be stars also are important in understanding the formation of massive stars. In 
additi on, it is possible that ac cretion plays a role in generating X-ray emission in HAEBE 
stars ( IHamidouche et al.ll2008l ). 



Hillenbrand et al.l (119921 ) suggested a classification scheme for HAEBE stars akin to 
that used for T Tauri stars. Using the slope of a star's spectral energy distribution in the 
near to mid-infrared, stars with large infrared excesses and spectral slopes with AF^ ~ A _4//3 
were assigned to Group I; these stars presumably have geometrically flat, optically thick 
circumstellar disks and may have optically thin inner disks. Group II stars have flat or rising 
spectra and were interpreted to have disk systems surrounded by gas and dust envelopes. 
Finally, Group III HAEBE stars have very small infrared excesses and appear to have only 
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gaseous emission from any residual circumstellar disk or envelope. 



The spectral energy distribution, which is composed of high spatial resolution optical 
and near-infrared observations combined with much lower spatial resolution Infrared Space 
Observatory (ISO), Infrared Astronomy Satellite (IRAS), and ground-based sub- millimeter 
continuum data may, however, be incorrect, if the far-infrared and sub-millimeter emis- 
sion seen toward these stars is dominated by emission from extended reflection nebulae or 
nea rby embedded protostars. I f this is the case, which it can b e, especially for B stars, see 
e.g.jHenning et al.l (119981 ll994j ): ISandell and Weintraubl (jl994al ); lAspin. Sandell k, Weintraub 



Jl994h and the results from this paper, then the disk and envelope models for these HAEBE 



stars must be re-examined based on higher resolution observations. 

In this paper, we present high spatial resolution 850 /im and 450 /im maps of HAEBE 
stars that enable us to examine what fraction of the long wavelength emission that has always 
been assumed to be part of the spectral energy distributions for these stars is associated with 
the stars and whether some of the emission originates from the surrounding cloud and/or 
invisible nearby young objects. 



2. Observations 



Almost all of the sub-millimeter observations reported in this paper have been obtained 
with SCUBA on JCMT0, Mauna Kea, Hawaii. Some of the observations were carried out 
by us in August 2001; others have been retrieved from the JCMT archive at the CADC0. 
We also present a few observations ob tained using the single element bolometer UKT14 on 
JCMT JSandell and Weintraubl fl~994bh . 



SCUBA ([Holland et al.lll999f ) has 37 bolometers in the long and 91 in the short wave- 
length array separated by approximately two beam widths in a hexagonal pattern. The field 
of view of both arrays is ~ 2'3. Both arrays can be used simultaneously by means of a dichroic 
beamsplitter. The long wavelength array is optimized for 850 /im and the short wavelength 
array for 450 /im, but one can also choose to observe at 750 /im and 350 /im for the long and 
short wavelength array, respectively, with slightly reduced sensitivity. The 750/350 /im filter 
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combination is no longer supported, but it was used in 1997 and 1998 and we have retrieved 
and reduced some 750 //m and 350 fim from this time period. The measured Half Power 
Beam Width (HPBW) is 14'.' 5 and 7'.'8 for 850 /im and 450 /im, respectively, for standard 
jiggle-map observing in stable night time conditions. For the 750/350 /im filter combination, 
we measured HPBWs of 12"9 for 750 /im and 6"5 for 350 /im, respectively. 

For sources with angular sizes less than or about equal to the field of view, the preferred 
observing mode is to use the long and short wavelength filter combination in jiggle-map 
mode. Larger emission regions have to be observed in scan-map mode. In jiggle map mode, 
the spacings between the bolometers are filled in by performing a jiggle pattern with the 
chopping secondary, while the secondary is chopping with a frequency of 7.8125 Hz. To 
completely sample both arrays, one needs to sample 64 positions with the array, while only 
16 jiggle positions are needed if only one of the arrays is used. The default time for each 
jiggle step is 1 sec for standard observing. Each 64-point jiggle pattern is done in a set 
of four, 16-point jiggles, each followed by a nod of the telescope into the other beam. In 
scan-map mode the array is tilted relative to the scan direction to ensure full sampling 
over the array with the telescope is scanning with a rate of 24" s _1 while chopping with 
a frequency of 7.8125 Hz. Each scan is overlapped by about half the arr ay. The t ypica l 



scan-map mode at JCMT is a basket weaving technique first suggested by lEmersonl (119951 ) 
where one can scan the source in an arbitrary angle, but chop in two orthogonal directions 
(fixed in equatorial coordinates) and restore the dual beam map in the Fourier space after 
converting the dual beam maps to equatorial coordinates. The standard setup for SCUBA 
is to use six maps, three of which are made while chopping in RA with three different chop 
throws and three while chopping in Dec. also with three different chop throws. The chop 
throws are chosen to so that the final map is sensitive to most spatial frequencies. Scan- 
maps never go as deep as jiggle-maps, because the restoration process results in baseline 
uncertainties and sky noise and eye n small pointing errors cause additional restoration noise 



(ISandell. Jessop and Jennessl 120011 ). However, for compact sources well above the noise level 
of the map, this is not a problem, since one can always subtract or filter out the background 
emission. 

The majority of the data presented in this paper have been obtained in jiggle-map 
mode, typically with a 120" chop in azimuth or at a fixed position angle in order to avoid 
chopping onto extended emission. A few compact isolated sources have been observed with 
shorter chop throws, but not less than 80". Most HAEBE stars are known to illuminate large 
reflection nebulae or are located in extended molecular cloud complexes; these sources have 
been observed in scan-map mode. A few stars — HD 200775 in NGC7023, and CoD— 42° 
11721 — have been observed in both scan-map and jiggle-map modes. For these regions both 
jiggle and scan-maps agree within errors with each other. The absolute pointing of SCUBA 
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is typically checked on blazars or calibrators near the target object before and after each map 
or set of maps. This enables us, in the data reduction stage, to correct for possible pointing 
drifts that may occur between the start and end of a map. Some of the maps retrieved from 
the archive did not have such frequent pointing observations. In these cases we have used 
the closest pointing observation to correct the pointing for the target star and then used the 
target star or another suitable point source in the map to correct for possible pointing drifts 
during the observations. The absolute pointing accuracy is expected to be better than 2". 
The sub-millimeter opacity was determined from sky dips at 850 /im and/or derived from 
CSO tau, the sky dip meter located at CSO, which provides opacity information every 10 
minutes at A = 1.3 mm. Two data sets were reduced from SCUBA in photometry mode. 
The reduction of of these data sets is essentially similar, although we only obtain a total flux 
density, i.e., no information about whether the source is compact or extended. 

All the data have been flux calibrated from maps of Mars, Uranus or Neptune or 
from maps of the SCUBA secon dary calibrator sources HLTau, OH231. 8+4.2, IRC+10216, 
RAFGL618 and RAFGL 2688 JSandelll fl994l Il998l : Ijenness et aDl2002h typically observed 
several times during a night. If sufficient calibration information was not available for a 
night, we have derived flux calibration from nearby nights with similar weather conditions. 
A complete list of the HAEBE stars observed in this survey is given in Table 1. 



The data were reduced in a standard way using SURF (IJenness and Lightfootl Il999 
Sandell. Jessop and Jennessl l200ll ) and STARLINK imaging software, i.e., we flat fielded, 
extinction corrected, sky subtracted, despiked, and calibrated the images in Jy beam -1 . We 
then pointing-corrected each scan for any drift in pointing between each pointing observation 
and added the data together to determine the most likely sub-millimeter position at 850 fim. 
Once we had derived a source position this way, we then applied small additional RA and Dec 
corrections to each scan (shift and add) to sharpen the final image to this position. Both de- 
spiking and sky subtraction can be done in several different ways depending on the amount of 
data and what observing mode was used for colle cting the data (ISandell. Jessop and Jenness 
200ll ) . We generally followed the advice given in ISandell. Jessop and Jennessl (120011 ) or tried 
several different ways in order to find the optimal way to ensure good data quality. 



All our maps were converted to FITS files and exported to MIRIAD (jSault. Teuben and Wright 



19951 ) for further analysis. In order to correct for the error lobe contribution, especially at 450 
/im, we have deconvolved all the maps using CLEAN and a circular model beam even though 
the beam is generally somewhat elliptical, i.e., broadened in the chop direct ion. However, for 
large data sets sky rotation will circularize the average beam profile anyway. Since our data 
set spans a variety of observing conditions and several observing seasons the beam size and 
error lobe contribution varies significantly from data set to data set. At 450 jttm the HPBW 
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can vary due to seeing conditions and changes in the surface by as much as 1". The error 
beam is even more sensitive to thermal gradients over the telescope. Especially in the early 
evening and during daytime, the 450 /iin error beam can change due to these thermal gradi- 
ents (Sandell, JCMT internal report). We find that the 450 //m HPBW varied by as much 
as 1" and occasionally more due to all these factors. 

For jigglemaps we measured HPBWs at 850 jttm in the range 14"3 - 15"5, but most of our 
maps in good stable night time conditions agree with the nominal value of 14'.' 5. At 450 fim 
the nominal HPBWis 7'.'8. For scan maps a typical HPBW at 850 /im is ~ 15", while it is ~ 
8" 5 at 450 /im. For jiggle-maps our model beam is composed of three symmetric Gaussians: 
a main beam, a near-error lobe, and an extended, low amplitude, far-errror lobe. For scan 
maps we only use a two component model, i.e., a main beam and a near-error beam. At 850 
jum the compact near-error beam has a HPBW of ~ 40" - 55" with an amplitude of 1 - 2% 
of the beam peak and an extended error lobe with a HPBW > 100" and an amplitude of ~ 
0.1% of the maximum. Because the amplitude of the extended error lobe is so small at 850 
fim, we can ignore it. At 450 /im we find the near-error beam to have a HPBW in the range 
25" — 34" with an amplitude 3 - 7% of the beam peak, i.e., much larger variations than 
those we see at 850 /im. The extended error lobe for jiggle- maps at 450 fim has a HPBW of 
~ 120" and an amplitude of < 1 %. 

In almost all cases, we have restored our CLEAN images back to 14" resolution at 850 
fim and to 8" resolution at 450 /im. In a few cases, where the signal-to-noise is marginal, 
we have restored the maps back to a resolution (15" - 16" for 850 /im; 10" - 12" for 450 
fim) which improves the signal-to-noise without degrading the resolution too much. Then, 
using the task IMFIT in MIRIAD, we fit a two-component elliptical Gaussian, one for the 
HAEBE star, and one for the surrounding cloud core. The fit of the cloud core is mainly 
used to provide a good subtraction of the extended emission and is not used to estimate the 
flux density of the surrounding cloud. For blended sources like LkHa 198 and LkHa 198 mm 
we used constrained fits with multiple components, increasing the uncertainty in the derived 
parameters. The results of these Gaussian fits for the detected HAEBE stars are given in 
Table l2al Note that the quoted flux densities are background corrected flux densities, i.e., 
the surrounding cloud emission has been subtracted before we computed the integrated flux 
density. Additional sub-millimeter photometric measurements with UKT14 are presented in 
Table [20 Limits for non-detections, all 3-cr upper limits, are given in Table [2cl 



-7- 



3. Archive Spitzer observations 



We searched the Spitzer archive for all MIPS 24 and 70 jj,m. observations of our target 
stars and retrieved and analyzed all photom etric (image) data. A few stars were covered by 
the c2d legacy project ( jEvans II et al.ll2003l ). For these stars we retrieved the results from the 



IRSA General Catalog Query Engine, Gator (http:/ /irsa.ipac. caltech.edu/applications/Gator/) 
If the MIPS pbcd-data appeared to be of good quality, we extracted photometry with Mopex, 
which was developed by the Spitzer Science center; otherwise we first reduced the images 
from the Archive BCD-data using Apex, another software package developed by the Spitzer 
Science center. The results of the Spitzer MIPS photometry are given in Table 



4. Analysis 



We fit the observed tota l fluxes with a simple graybody fit, see e.g., ISandelll (120001 ); 
Sandell and Weintraubl (120011 ) . For "isolated" HAEBE stars, i.e., stars for which we do not 
see any emission from the surrounding cloud or from nearby sources, we make use of IRAS, 
ISO or KAO far infrared data longward of 60 /im in addition to our own sub-millimeter data. 
In a few cases, where we have no 100 /im data or where there is a clear excess at 100 /im, 
we assume a dust temperature or explore a plausible range in dust temperature. The fits for 
these stars are shown in Figure [1] and Figure [2j 

For confused regions, i.e., where nearby sources or strong emission from the surrounding 
cloud core are present, we do not use the far-infrared data, because we have no way to 
accurately partition the contribution of each component to the far infrared flux measured 
in a much larger beam. In these cases we also assume a plausible dust temperature based 
on all the available information we have for for these sources.. Since we have measured the 
source size (Table I2"a|) . we constrain the size and only fit for the dust emissivity, (3, the dust 
temperature, T^ (if we have far infrared flux density measurements), and the dust optical 
depth at 850 /im. To derive masses we as sume standard Hi ldebrand opacities (ki2oogh z = 0.1 



cm 2 g x ) and a gas-to-dust ratio of 100 (IHildebrandl 1 1 9831 ) . The results of the graybody fits 



are given in Table HI In the same table we also list the spectral index of the dust emission 
a (F u oc v a \ after subtraction of free-free emission, which was derived by least-squares 
fitting all millimeter and sub-millimeter data for each source. For optically thin dust in the 
Rayleigh- Jeans regime, the spectral index, a = (3 + 2, where j3 is the dust emissivity index. 
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5. Results 

All maps of the sources listed in Table 1, except one, show emission in the SCUBA 
field of view, although the observed emission is not always associated with the HAEBE 
stars. Early B stars are luminous enough to heat up their surrounding clouds and, with the 
exception of the very distant HAEBE star V431 Set (MWC 300), invariably show extended 
cloud emission. This is the only field where we see no emission at all. We have only detected 
three HAEBE stars of spectral type early B or late O. The failure to detect early type 
stars does not, however, mean that massive stars do not form disks; rather, in this sample 
it appears to be a selection effect. When Herbig Be stars are bright enough to illuminate 
reflection nebulae, they may already have freed themselves from the cloud cores in which 
they were formed, i.e., they are already close to the main sequence and appear to have 
dispersed the disks with which they were born. 



5.1. Isolated HAEBE stars 



Isolated HAEBE stars have the same characteristics as "classical" HAEBE stars, i.e., 
they are often located near, but not in star forming regions, they have Ha in emission, IR ex- 
cesses due to thermal emission from dust, but differ from them due to the absence of nebulosi- 



ties (jGrinin et al.lll99ll ). The absence of nebulosities suggests that these stars on the average 
are older than "classical" HAEBE stars, because they have already freed themselves from the 
clouds in which they were born. At least three (HD 135344 B, H P 141569, and HP 169142) 
possi bly more of our isolated HAEBE stars are transition objects (INajita. Strom &: Muzerolle 
20071 ). These are stars which show significant disk evolution. The transition disks are char- 
acterized by an inner gap, which is still gas rich but contains very little dust, and a cold 
outer disk. Pue to the inner gap they have very little near-infrared excess, but they have 
strong mid- and far-infrared excess due to the outer disk. 



Meeus et al.l (120011 ) did ISO spectroscopy of 14 isolated HAEBE stars and separated 



them into two groups based on their infrared SEPs: group I and group II. For group I the 
infrared to sub-millimeter continuum could be modeled by a power-law and a blackbody, sug- 
gesting that stars of this group have optically thick and geometrically thin disks (power-law 
component) and an optically thin flared region (black body compon e nt). G roup II objects 
only needed a power-law to fit their continuum, which iMeeus et al.l (120011 ) suggested indi- 
cates that they are more evolved stars with partially optically thin inner disks protecting 
the outer disk from flaring. We have sub-millimeter data on many of the stars observed by 
Meeus et al.l (120011 ) . Observations of several of them made with millimeter arrays confirm 
that the infrared excess is due to circumstellar disks ([Mannings. Koerner. fc Sargentlll997l ; 
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Mannings fc Sargentlfl997l . I2OO0I : iNatta et al.lbooi lHamidouche et al.lbood : llsella et al.f2007h . 
In the millim eter and sub-millimeter regime there is very little difference, if any, between 
Meeus et al.l group I and group II objects. 



tral class of A3 - A5 Ve and an age of 4 



(van den Ancker et al. 


1998). It has a spec- 


8 Mvr ( 


Simon, Dutrev & Guilloteau 


2000; 



Pietu et al.l 12006c iPietu. Dutrey fc Guilloteau! 120071 ). The star has never been searched for 
free-free emission, but if it has free-free emission, it is presumably weak. It has been exten- 
sively stu died with mm-arrays and is surrounded by a large circumstella r disk in Keplerian 



rotation (Mannings fc Sargent 



Pietu. Dutrey fc Guilloteau 



1997 



2007 



Simon. Dutrev fc Guilloteau! I200CH : IPietu et al.l 12006 



Hug hes et al.l l2008af). The disk is well resolved both in 



CO and continuum. At 1.4 mm lHamidouche et al.l ( 120061 ) find a disk size of 0''8 x 0''7 at a 
P.A. of 143° ± 5°, while the CO disk is much larger with a diameter o f ~ 8" at a P.A. of 



150° rtSimon. Dutrey fc Guilloteau! |2000|; IPietu. Dutrey fc Guilloteau! 12007). Both contin 



uum and CO images give a similar incli nation angle for the disk, i = 37° (lHamidouche et al. 



20061 ; IPietu. Dutrey fc Guilloteau! 120071 ) . Even though the di sk is quite large, it is not seen 
i n cor onagraphic near-infrared images ( Augereau et al.l IboOll ) . IPietu. Dutrey fc Guilloteau 
(120071 ) derive ~1.8 M for the mass of the central star from analysis of the Keplerian rota- 
tion of the disk in the J = 1 — > and J = 2 — >■ 1 transitions of CO and its isotopomers as 
well as from HCO + J = 1 — )> 0. The derived mass would suggest that the star is of spectral 
type A5. 

Our SCUBA observations (Table l2al Figure [3]) resolve the disk around MWC480 both 
at 850 and 450 fim. Our observations show an emission feature extending to the south 
from the disk. Thi s emis sion is nearly at the same P.A. as the observed jet-like emission by 
Hamidouche et al.l (120061 ) and could therefore be real. Our isothermal graybody fit to the 
millimeter and 100 fim data predict a cold dust disk, T^ ~ 28 K and a dust emissivity index, 
(3 ~ 0.8 (Table H]), As we can see from Figure [1] the dust emissivity is well constrained, but 
the fit underestimates the observed flux densities at 60 fim, which suggests that one needs a 
warm dust component as well in order to explain the 60 /im emission. 



HD 34282 is of spectral type A3 Ve flMora et al 



2001 



Merin et al.l 12004! ) with an age 



of 6.4 Myr. It has no detectable free- free emission (jNatta et al.l 120041 ) . Its disk was resolved 
at 1.3 mm with the IRAM Pla teau de Bure interferometer (PdB) with a size of 1"7 x 0"9 
(IPietu. Dutrey fc Kahandl2003l ). They find from imaging the disk in CO J = 2 — > 1 that the 
disk is in Keplerian rotation and seen at an inclination of 56° ± 3°. Based on the dynamical 
mass they argue that the previously reported distanc e estimate, from Hipparcos parallax 
observations, 160 ± 42 pc (Ivan den Ancker et al.lll998l ) must be in error and that the most 
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l ikely distance is 400 pc. The detailed modeling of the stellar parameters by iMerin et al. 
( 120041 ) gives 348 pc, We will ther efore use the distance 350 p c. The disk was also imaged 
with th e OVRO interferom eter by lMannings fc Sargentl (120001 ) in the continuum at 2.6 mm 
and by iNatta et al.l (J2004J) with the PbB inter ferometer at 3.2 and 1.3 mm. T he observed 
continuum fluxes agree well with the results of iPietu. Dutrey fc Kahand (120031 ) . 



We find the star to be un resolved at 850 /im and 4 50 um (Table [2a|) . These SCUBA data 
have also been published by ISheret. Dent fc Wyattl (120041 ). who get somewhat discrepant 
results at 450 /zm. In this analysis we use our result s. The star was included in the CO J = 
3 — )► 2 single dish JCMT survey by lDent et al.l (120051 ) and shows a CO profile similar to what 
was observed with the PdB interferometer, confirming that this is an isolated HAEBE star. 



Since the star appears compact in both single dish, coronographic imaging (IDoering et al. 



20071 ) and aperture synthesis observations, we combined all millimeter and sub-millimeter 
data in our isothermal modeling. The d ust emissivity we derive, (3 = 1.29, is very similar 
to what IPietu. Dutrey fc Kahand (120031 ) derived from their detailed modeling. However, in 
this case we find the dust emission to be marginally optical ly thick at 850 fim, r^n ~ 0.2, 
and we get a slightly higher mass, 0.18 M Q (Table H]), than IPietu. Dutrey fc Kahand . The 
mass of the disk is much higher than we find for any other isolated Ae star. 

HD 35187 is a visual binary s tar with a separation of 1'.'38 in Taurus at a distance of 
150 pc (IDunkin fc Crawfordlll998l ). iDunkin fc Crawford! found the stars have spectral types 
of A2 V (HD 35187 B) and A7 V (HD 35187 A). Only HD 35187 B shows Ha in emission and 
has an infrared excess, suggesting it is surrounded by an accretion disk. The age of the 
system is ~ 10 Myr. 



HD 35187 was detected at 1.1 mm and 800 fim by lSylvester et al.l (119961 ) with the single 



(Natta et al. 


2004 


)• 


Natta et al. 


(2004) 



with the VLA at 7 mm and with the PdB interferometer at 3.6 and 1.3 mm. Most of the 
emission at 7 mm is due to free-free emission. The free-free emission is still significant at 
3.6 mm, but at 1.3 mm and shorter wavelengths the emission is dominated by dust emission. 

Our 850 and 450 /im SCUBA maps show an unresolv e d poi nt source. These SCUBA 
maps were previously published by ISheret. Dent fc Wyattl (120041 ) . We have reduced them 
independently and find somewhat higher flux densities (16% and 27% at 850 fim and 450 
//m, respectively) than ISheret. Dent fc Wyattl . but in good agreement with other published 
data on this star. A graybody isothermal fit to all (sub-)millimeter data and the IRAS 100 
/im data point (Figure [Q gives a dust temperature of 44 K and a dust emissivity index, j3 
= 0.87. The mass of the disk is relatively low, 0.002 M (Table 0}. 
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HD 3 6112 is an isola t ed HA EBE star located in the outskirts of the Taurus-Auriga 
complex. IChapillon et al.l (120081 ) adopted a distance of 140 pc, arguing that the star is 
likely to be at the same distance as the Taurus- Auriga complex. iBeskrovnaya et al.l (119991 ) 
classified it as A8 Ve with an age of 5 - 10 Myrs. A three component fit to low resolution 



UV spectra give a spectra 



pc ([B londel & Tiin A Diie 



(Ivan den Ancker et al 



class of A7 Hie, with an age of 4.4 Myr and a distance of 200 



20061 ). The latter agrees with the Hipparcos parallax distance 



19981 ). We assume 200 pc. 



The star ha s not been detected in free- free emission at centimeter wavelengths. It was 
first observed by iMannings &: Sargentl (119971 ) with the OVRO array in continuum at 2.7 mm 
and in the 13 CO J = 1 — >• line. The dust disk was unresolve d in continuum and barely 
detected in 13 CO. Follow-up observations with the OVRO array (IMannings Sargentl 2000) 



observed the star in continuum at 1.3 mm, but did not resolve it. IChapillon et al.l (120081 ) 
observed the star with the IRAM PdB interferometer in continuum at 3.4 mm, 2.6 mm, 
and 1.3 mm, as well is the 12 CO J = 2 — > 1 and 1 — > transitions. They resolve the disk 
in continuum at 1.3 and 2.6 mm with a size of ~ 1" (1 50 AU), and find the disk slightly 
more extended in CO (r ~ 250 AU). bent et al.l (boosh detected the disk in the CO J = 



3 — » 2 line emission. By modeling the CO profile they found the disk to be almost face-on 
with an inclination < 10° and deduced an outer radi us of 170± 30 A U. Their results agree 
quite well with the far better constrained modeling by IChapillon et al.l . who determine a disk 
inclination of 16° - 18° from their CO observations. 

We do not resolve the dust emission, which agrees well with the the PdB results. A 
graybody fit to all (sub-)millimeter and FIR data to 100 /iin, including ISO data from 
Elia et al.l (120051 ) . give a dust emissivity index, (3 = 1.25, a dust temperature of 43 K, and a 
disk mass of 0.014 M Q . Our deduced dust emissivity is slightly higher, but well constrained 
by our fit (Figu r e [1]), and our disk mass is a factor of two lower than what was found by 
Chapillon et al.l ( 120081 ) . They derived a dust emissivity index, (5 = 1.0 ± 0.15, and a disk 
mass of 0.027 M . 



HD 135344 B is an isolated F4 Ve star ( IDunkin. Barlow, fc Ryan 



2009) in the Sco OB 2 - 3 association at a distance of ~ 140 pc (ICoulson fc Walther 



1997 



Grady et al. 



1995 



van Boekel et al.l 120051). The star has been extensively studied from X-rays to millime- 



ter wavelengths (iGradv et al.l 12009k and references therein). It has strong millimeter and 



far-in frare d excesses dWalker fc Wolstencroftlll988t ICoulson fc Waltherlll995t ISylvester et al. 



19961 ). and lGrady et al.l (120091 ) co nclude that is definitely a PMS star w ith an age of ~ 8 Myr. 



As deduced both from modeling (IDent et a 



the mid-infrared and in the sub-millimeter ( IPontoppidan et al. 



2005; 



Brown et a 



200 7) and direct imagi ng in 



20081 : iBrown et al.ll2009l ). the 



disk is seen approximately face-on and has a large inner gap. The star therefore has a tran- 
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sitional disk, where the in ner disk is largely free of dust. Yet the disk appears surprisingly 
gas rich (IDent et al.ll2005l ). 



The st ar is well studied in the (sub-)millimeter wavelength regime. Both lCoulson &: Walther 

rli995h and ISvfvester et al.l (119961 ) obtained photometry of the star with UKT14 on JCMT, 



Walker fc Butnerl (11995T) observed i t at C altech sub- millimeter Observatory (CSO) at 1.3 



mm and 800 /mi, and iBrown et al.l (120091 ) imaged it with high spat ial resolution with the 
Smithsonian Submillimeter Array (SMA) at 880 jttm. IBrown et al.l resolved the dust disk 
with a radius of ~ 0'.'9 (125 AU) and found an inner hole with an outer radius of 39 AU. 
We have deep images with SCUBA at 850 /mi and 450 /im, which show that the disk is 
unresolv ed in an 8" beam (Table l2a|) . This is consistent with the observed size (1'.'8) at 
880 /im (IBrown et al.l 120091 ) . although our observed flux density is ~ 50% higher than what 



Brown et al.l find in their high spatial resolution SMA observations, suggesting they may 



filter out some of the extended emission. Our isothermal graybody fit, including all mil- 
limeter, sub-millimeter as well as the IRAS 100 /im flux densities gives a dust emissivity 
index, j3 = 1 .38, and a dust temperatur e of 32 K. We did not include the ISO 200 /mi 
flux density (jWalker fc Heinrichsenl 120001 ) . because it appears anomalously low compared 
to all the sub-millimeter a nd IRAS data. We derive a much lower dust temperature than 



Coulson fc Waltherl (119951 ) . because they fitted the SED with a blackbody. As with all 



evolved disks, however, the graybody fit to the millimeter and sub-millimeter data under- 
estimate the flux densities shortward of 100 /mi, which require a warm dust component as 
well. If we assume normal gas-to-dust ratio, we get a disk mass of 0.025 M (Table HJ). 

HD 141569 is a nearby, isolated AO Ve/B9.5 Ve star, which is the primary mem- 
ber of a triple system in the outskirts of the Sco-Cen s ystem with an age of 3 - 5 Myrs 
(IDunkin. Barlow, fc Ryanl 119971 ; IWeinberger et al.l l2000h . iMerin et al.l (120041 ) determined a 
more precise age for the primary by fitting surface gravity, metallicity and effective tempera- 
ture to high resolution spectra and photometry. They find an age of 4.7 Myr, a spectral type 
of B9.5 Ve and a distance of 108 pc, which agrees wit hin errors with the Hipparchos distance, 



99 pc . Coronagraphic imaging in the near- infrared (jAugereau et al.lll999t IWeinberger et al. 
19991 ) shows a large (r ~ 400 AU ), slightly asymm etric dust disk with depression (gap) at 



250 AU and a 30 AU inner hole (Marsh et al. 



mid-infrared (IFisher et al.l |2000| ; iMarsh et al. 



2002). The disk has also been resolved in the 
faooj) with a size of 2^'2 at 10 /mi. HD 141569 
is often referred to as a debris disk; however, IDent et al.l (120051 ) detected CO J = 3 — > 2, and 
all debris disks are gas poor. The only su pposed debris disk star known to have a molecular 
gas disk is 49 Ceti (IHughes et al.ll2008bl ). and it is quite possible that the disk surrounding 
49 Ceti is primordial, rather than a debri s disk. In addition Spitzer IRS spec tra also show 
that HD 141569 has strong PAH emission (ISloan et alihood ; Keller et al.ll2008l ). which is not 
seen in debris disks. Furthermore, at an age of only 5 Myr, HD 141569 is likely too young to 
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be a debris disk star. Thus, its youth, the presence of a gas in the disk, and the detection of 
PAH emission from the disk all strongly suggest that HD 141569 has a primordial transition 
disk, not a debris disk. 

For this star the the SCUBA data presented in Table l2alwere obtained in photometr y 



mode and are based on the same data that were published bv lSheret. Dent fc Wvatt 



The flux densities that we obtain agree quite wel l with those of ISheret. Dent fc Wvattl . The 
star was also observed by SCUB A at 1.35 mm by ISvlvester. Dunkin fc Barlowi (120011 ) and at 
CSO and IRAM at 1.3 mm by (IWalker fc Butnerlll995l ). but the latter have poor signal-to- 
noise and have not been used in our analysis. The star was observed in photometry mode 
by Spitzer MIPS at 24, 70, and 160 fim. The photometry at 24 and 70 /xm is given in 
Table Eland agrees with results by Weinberger (private communication). The MIPS images 
only show emission from HD 141569 and no emission is seen from the companions. Our 
graybody fit severely underestimates the observed flux at 100 and 70 fim, suggesting that 
the disk must have a warm dust component as well. The temperature of the cold outer disk 
is therefore uncertain. In Table H] we give the results for an assumed dust temperature of 40 
K, which results in a f3 index of 0.5 and a disk mass of 9xl0~ 5 M . Because we assumed 
a temperature rather than derived it, this results in much higher uncertainties for the dust 
emissivity and disk mass. By varying the temperature from 25 - 70 K, we find that varies 
by less than ± 0.2, while the mass could be off by a factor of two. 



HD 142666 is located in the Sco OB2 -2 association at a distance of 145 pc (Ivan Bqekel et al. 



20051). It has a spectral type of A8 Ve ( Dunkin. Barlow. &: Ryan 1997 ; Mora et al. 



2001 



van Boekel et al. 2005 ; Guimaraes et al. 20061 ) with an estima ted age > 10 Myr ( Natta et al. 



2004 ). It was not detected in free-free emissio n at 3.6 cm ( iNatta et al.l I2004J ) . It was in- 
cluded in the ISO survey by lMeeus et al.l ( 120011 ) where it was classified as a group Ha object, 
i.e. a star having an optically thick disk in the FIR. The disk was detected in the CO J = 
3—^2 survey by ( iDent et al.l 120051). E ven though the CO emission is clearly detected, the 
signal-to-noise is poor and IDent et al.l fitted it with a single broad line predicting that the 
disk is seen nearly face-on (i = 18 ± 5°) and a radius of 45 AU, which would suggest a disk 
size < 1". 



The star was observed with a 7-channel bolometer array on IRAM at 1.2 mm (iBockelee-Morvan et al 



1994) and with UKT14 on JCMT at wavelengths from 2 mm to 450 /mi (Ivan der Veen et al. 



19941 ; ISvlvester et al.lll996l ). These observations were all photometric observations. HD 142666 
was also observed with the IRAM PdB interferometer in continuum at 3.3, 3.1, and 1.2 mm. 



Although INatta et al.l (120041 ) give no size information, the 1.2 mm flux density from PbB 
appears somewhat low compared to single dish photometry, suggesting that the disk could 
be extended. The disk is unresolved in our SCUBA images at 850 /xm, although at 450 fim 
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we obtain a size of 2'/2. Since the size is not verified at 850 /jm, we consider this an upper 
limit. Far infrared flux densities observed with ISO (jElia et al.ll2005l ) appear overestimated 
when compared to IRAS data, and we have not included them in our analysis. If we only use 
the (sub-)millimeter observations and the IRAS 100 /im flux density, our graybody fit gives 
a relatively low dust emissivity, f3 = 0.72 for a dust temperature of 31 K. The isothermal 
fit cannot reproduce flux density at 60 /im (Figured]). To fit the millimeter and far- infrared 
SED we need a warm dust component as well, which is true for most isolated HAEBEs. 

HD 144432 is a binary star w ith a late type c ompanion in the Sco OB2-2 asso- 
ciation at a distance of ~ 145 pc (IPerez et al.1 [200J). The spectral type i s somewhat 
uncertain. iDunkin. Barlow, fc Ryan! (19971 ) classified it as A9/F0 Ve, while IPerez et al. 



J2004h and blondel fc Tiin A Due! (boog T assign it a spectral stype A8 Ve or A9 Hie with 



an age of 6.5 - 6. 8 Myr. It was included in the ISO survey of isolated HAEBE stars 
(IMeeus et al.1 120011). where it was c lassified as group Ha object. It was detected at 1. 2 
mm by iBockelee-Morvan et al.1 (119941 ) and at 1.3 mm and 800 /im by lSylvester et al.l (119961 ). 



We have only SCUBA photometry of this star (Table I2"a|) . A graybody fit to all (sub- 
millimeter and IRAS 100 //m data give f3 = 0.59 and a disk mass of 0.002 M Q . The fit 
greatly underestimates the flux density at 60 /im (Figure [2]), suggesting that we need a warm 
dust component as well. 



HD 150193 is located in Op hiuchus at a distance of ~ 1 50 pc. iMora et al.l (120011 ) 



a spectral type F9Ve. 
the binary system of r 



classified it as an A2 IVe star, while iBlondel fc Tiin A Diid (120061 ) found it to have a spectral 
type A3 Hie from analysis of UV-spectra. Eliaa (1978) cla ssified it as AO. It is a visual 
binary with a separation of I'll (jReipurth fc Zinneckerl 1993 ). The companion is quite faint. 
Carmona. van der Ancker fc Hennina (120071) find the secondary to be a T Tauri star with 



Carmona. van der Ancker fc Henningi (120071 ) find a common age for 
10 Myr. The star was not detected in the CO J = 3 — » 2 survey by 
Dent et al.l (120051 ). confirming that it is an isolated HAEBE star, which is not associated with 
any molecular cloud. It also appears to have a gas-poor disk. HD 150193 is associated with 



faint free-free emission, 0.2 mJy at 3.6 cm (ISkinner. Brown, fc Stewartl Il993l ; iNatta et al. 
20041 ). almost certainly originating from a thermal wind. 



The star has strong infrared excess (lEhaslll978l ) and was first detected in sub - millimeter 



contin uum at 800 and 450 /iin with the single pixel bolometer UKT14 on JCMT by lJensen. Mathieu fc Fulle 
(119961 ) . It was imaged wi th the OVRO interferomete r in the 2.6 mm dust continuum and in 
the CO J = 1 — > line ([Mannings fc Sargentlll997l ). who found t he dust continuum to be 
unresolved and did not detect the disk in CO. It was also imaged by lNatta et al.l (120041 ) with 
the VLA at 7 mm in their study of evolved dust disks in HAEBE stars. We report SCUBA 
observations at 850 /im (Table l2aj) . which show an unresolved point-source. The star was 
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also imaged by Spitzer MIPS, which show an isolated point source at 70 /mi (Table [3]). It 
was not detected by IRAS at 100 /mi, but at 60 /mi the IRAS flux density is more than twice 
as large as what we see with MIPS at 70 /mi suggesting that IRAS picked up some extended 
emission in its much larger beam. Fitting an isothermal graybody disk model to all available 
millimeter and sub-millimeter data and the MIPS 70 /im data point indicates that the disk 
is marginally optically thick at 850 /im (r ~ 0.2) with a temperature of 44 K (T able HI) and a 
mass o f 0.002 M Q . The dust emissivity we derive, (3 = 0.65, is lower than what iNatta et al. 
(120041 ) found ( j3 = 1.0 - 1.6) from their more limited sampling of the SED. Figure [2] shows 
that our fit is reasonably well constrained with an uncertainty in j3 of ± 0.2. 

KK Qph is in the o utskirt s of the Ophiuchi cloud comp lex. It was classified as A8 
Vev bv iMora et aD (boOlh . while iBlondel fc Tim A Duel (120061 ) classify it as a A5 Ve, with 
a range A5 - A7 Ve from modeling UV spectra and photometry. It is a binary with 
a separation of F /5 (IPirzkal et al.l 1 1997) The secondary is a T Tauri star with a spec- 
tral type G6 Ve (ICarmona. van der Ancker fc Henningi 120071 ). The distance is not we ll 
deter mined, but IBlondel fe Tiin A Diie J2006h quote < 250 pc from IWallenquistl (119371 ). 



while Leinert et al. 



(20 04]) estimate the distance to be 160 ± 30 pc, which we adopted in 



this paper. IBlondel fc Tiin A Pi id es timate an age of ~ 7.9 Myr, which agrees well with 
Carmona. van der Ancker fc Henningi (120071 ) who find an age of 7 Myr. Even though Spitzer 
MIPS archive images at 24 and 70 /im do not show any sign of nebulosity or cloud emission 
around the star, the flux density observed by MIPS at 70 /im is rather faint when compared 
to the IRAS 60 /im flux density. The MIPS flux density at 70 /mi, even without color correc- 
tio n is only 2.88 ± 0.02 Jy , while the ADDSCAN IRAS flux density at 60 /im is 6.7 ± 0.7 



Jy ([Hillenbrand et al.lll992l ). Since both IRAS and Spitzer produce well calibrated data, it 
would appear that the much broader IRAS beam picks up some e xtended emission around the 



star. The star has weak free-free emission, 0.25 mJy at 3.6 cm ( jSkinner. Brown, fc Stewart 



1993|). 



KK Oph was detected 1.3 mm by ( [Hillenbrand et al.lll992l ) with broad beam photom- 
etry, but there are no other published (sub-)millimeter observations. We obtained 800 /tm 
photometry with UKT14 (Table l2bl) and mapped it with SCUBA at 850 /mi. We find 
KK Oph to be faint and unresolved (Table I2"aj) . For KK Oph, we derive a fairly typical mass 
of 0.002 M and a rather shallow f3 index, 0.74; however, the j3 index is not well constrained, 
since we only have data covering a narrow wavelength regime. 

HD 163296 is an early A star at a distance of 122 pc and with an age of ~ 4 Myr 
(Table [T]). The star is isolate d, i.e., there is n o sign of th e molecular c l oud in which the 
star was presumably formed ( IDent et al.ll2005l ). However, iGrady et al.l (120001 ) show that 
HD 163296 is still an active Ae star, suggesting that there are still some remnants left of the 
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cloud out of which it was formed. Single-dish observati ons of HP 163296 show strong dust 



emiss i on at millimeter and sub- millimeter wavelengths ([Hillenbrand et al.lll992l; Mannings 



1994 ; iHenning et al.l Il994j . Il998l ). The star is very faint in free- free emission ( iNatta et al 



20041 ) . Therefore, the m illimeter emission is completely dominated by thermal emission from 



dust. Ilsella et al.l (120071 ) have observed HD 163296 in continuum with the SMA at 870 /im, 
with PdB at 1.3 mm and 2.8 mm, and with the VLA at 7 mm. They resolve the disk at 
all wavelengths except 7 mm, where the VLA is less sensitive or has resolved out the cold, 
extended, outer part of the disk. At 870 /im they derive a disk size of 3''6 x 2"6. They also 
imaged the disk in emission from several molecular lines. From the 12 CO J = 3 — > 2 and 
12 CO J = 2 — >• 1 lines, they deduce that the velocity pattern is well fitted with a Keplerian 



disk w ith an inclination angle 46° ± 4° and a stellar mass of 2.6 M , while iGrady et al. 
( 12 00 Oh found an inclination angle of ~ 60° ± 5° from their coronagraphic imaging of the 
reflected light from the disk. The radius of the disk observed in CO, ~ 540 AU, is much 
larger than in continuum, ~ 200 AU; such a behavior can be explained by an exponential 
falloff of the surface density at large radii ( jHughes et al.ll2008ai ). 

We find HD 163296 to be extended and well resolved at both 850 and 450 /im (Table 
I2"a] Figure H]). Although there are many observations of HD 163296 in the JCMT archive, 
very few of the observations were done in good sky conditions with sufficient calibration 
observations to enable us to obtain well-calibrated data sets with good beam characterization. 
We retained only three data sets taken in engineering time in 2004 and 2005. Even though 
we find the dust emission around HD 163296 to be slightly elliptical, the elongation is not 
large enough to permit us to determine the position angle of the disk with high accuracy at 
either 850 or 450 /im. For our 450 /im data sets, where the HPBW is ~ 8", we created a 
model beam by rotating the measured beam profiles to the average position angle of each 
map of HD 163296, since the beam is always somewhat elongated in the chop direction 
(see Section [2]). We added the three rotated beam profiles together to form a beam profile 
which best approximates the beam of the coadded 450 /im map. After cleaning the 450 
/im map with this model beam, we find a FWHM = 4"6 x 3"2, P.A. = -13° ± 27°, with 
an u ncertainty o f ~ 1" (Table l2a|) . The size of the 450 fim disk is slightly larger than 



what Ilsella et al.l (120071 ) determined from their 870 /im SMA observations, but both re sults 
agree well within observational errors. The total flux densities derived by Ilsella et al.l also 
agree very well with single dish results at 1.3 mm and 870 /im, confirming that the aperture 
synthesis observation have recovered all the continuum emission from the disk. We can 
therefore use both interferometric data and single dish data to determine the emissivity 
index, f3, and disk mass. An isothermal graybody fit to these data gives T = 28 K and j3 
= 0.94. By assuming a gas-to-dust mass ratio of 100, we deduce a disk mass ~ 0.065 M Q . 
The graybody fit underestimates the IRAS flux density at 60 /im, suggesting that there is 
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an additional contribution from warm dust at 60 /an, although its cont ribution to the dis k 
mass is negligible. Our results are very consistent with the results of llsella et all ( 120071 ). 
These studies indicate that some grain growth has already occurred in the HD 163296 disk. 



HD 169142 was classified as A5 Ve (Dunkin. Barlow, fc Ryan 1997; van Boekel et al. 



2005 ). A7Ve falondel fc Tim A Diiei[2006h and as A9 IH/IVe JGuimaraes et al-bood Ivieira et~al 
20031 ). Here we adopt the spectral class A7 Ve, which is based on fitting low resolution 
UV spectra with a three component model, the stellar photosphere, an optically thick 
accretion disk, and a bound ary layer between the disk and the star. The analysis of 
Blondel fc Tjin A Djid (120061 ) indicate that the disk is seen approximately face-on, which 
agrees with agrees with other findings (see below) and predicts a distance of 150 pc and 
an age of ~ 12 Myr . Thi s distance agrees very well wit h the distance of 145 pc derived 



by Ivan Boekel et al.l (120051 ). which we have adopted here. iDent et all ( 120061 ) did not detect 
any free-free emission at 3.6 cm, although the flux density they observed at 7 mm appears 
anomalously high and is most likely due to free-free e mission or a radio flare (see below). It 
was included in the ISO survey of isolated HAEBEs (IMeeus et al.l 120011) and classified as a 
lb type HAEBE s tar. HD 169142 sh ows a narrow emission line in CO (IDent et al.ll2005l ) and 
in [Oil at 6300 A(Acke et al. 1 120051 ) . suggesting that the disk must be seen nearly face-on. 
Meeus et al.l (120101 ) modeled the disk including new Herschel data, and found that in order 
to explain the weak 10 /im excess emission, the disk must have a gap, i.e. it appears to be 
a transition disk. 

In the (sub-)millimeter the dis k was first detected in photom etry mode at 1.2 mm with 



the 7-channel bolometer on IRAM 



Bockelee-Morvan et al.lll994l ) and at 2 mm, 1.3 mm, and 



800 /im with UKT14 on JCMT ([Sylvester et alll996l ). Our observed continuum flux densities 
(Table l2aj) agree well with the phot ometry results and s uggest that the disk may be resolved 
with a size of 2'.'3 x 1"6 (Figure E). iRaman et aD J2006h observed HD 169142 with the SMA 
in continuum at 1.3 mm and in the CO J = 2 — > 1 line. They barely resolve the continuum 
emission (FWHM ~ 0'.'9) but find that the CO emission traces a Keplerian disk around the 
star. Their best fit disk model gives a radius ~ 235 AU (1''6), a disk mas s of ~0.02Mp, an d 
an inclination of 13°. This result agrees well with the SED modeling by lDent et al.l ( 120061 ). 
who find a disk radius of 300 AU, a disk mass of 0.04 M and an inclination angle of ~ 
30°. Our observed disk size, r ~ 140 AU, is consistent with these previous measurements. 
The star was also obse rved by the PACS imager on Herschel in continuum at 160 and 70 
/xm (IMeeus et al.ll2010l ). Our graybody fit inclu des all (sub-) millimeter and far infrared data 
up to 100 /im, except the 200 /im ISO value ( IWalker &: Heinrichsenl 120001 ) . which appears 
anomalously low, and the 7 mm VLA data point, which is too high to be due only to dust. 
It must therefore include a significant contribution from free-free emission. When we include 
the accurate 160 /im flux density from PACS, we get (3 = 1.30, a dust temperature of 41 K, 
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and a disk m ass of 0.03 (Table H] and Figure [2]). The latter is in good agreement with 



the results of iRaman et all. The dust em issivity that we obtained, j3 = 1.30, is more than 



twice as large as what iDent et all (120061 ) found, because they did a fit to a more limited 
set of data and assumed that all the flux density at 7 mm was due to dust emission. As 
we can see from Figure [21 the fit underestimates the observed flux densities shortward of a 
100 jum, requiring a warm dust component as well. Since we have reliable far-infrared data 
on HD 169142 we also fitted the the millimeter and far-infrared SED with a two-component 
graybody model, where we allowed the dust properties, dust temperature and size of the 
emitting region to vary for each component. However, we set the dust emissivity for the 
warm component equal to two, although the fit is not very sensitive to the emissivity of the 
warm gas. The two-component fit results in a slightly higher dust temperature for the cold 
dust, 52 K vs. 41 K. The dust emissivity for the cold dust component is now 1.38 compared 
to 1.30 for the single temperature fit, i.e., a marginal change. The temperature of the warm 
dust is 155 K. The two-component graybody fit now matches the 70, 60, and 25 /im flux 
densities very well, see Figure [2j 

MWC 349 is a remarkab le binary system at a distance of 1.2 kpc wit h a luminosity of ~ 
3 x 10 4 L JCohen et al.lll985h . MWC 349 A was classified as a Be star bv lMerrill fc Burwell 
(119331 ) . Its opti cal and infrared spectrum is complicated. T he star has very strong hydrogen 



and He I lines ([Thompson et al 



19771 ; I Cohen et al.1 119851 ) . Many lines are double peaked 



sugge sting they are formed in a Keplerian disk viewed edge-on (IHamann fc Simon 



1986 



1988). The secondary is at a separation of ~2V4 and has a spectral type B0 III (ICohen et al. 
19851 ). MWC 349 A is the brightest known radio star at centimeter wavelengths and the 
only known hydrogen recombination line maser source. Maser action is seen in hydrogen 
recombination lines from ~ 100 GHz (H39a) to 1800 GHz (H10a), with the all masering 
lines being double peaked and with the two p eaks coming from the receding and app roaching 
parts of the disk, due to Keplerian rotation (IStrelnitski. Smith fc Ponomarevl 119961 ). Direct 
imaging of the H30a recombination line with the SMA at 231.9 GHz confirm that the maser 
spots originate in an approximately edge-on disk with a r otation curve that is s omewhat 
steeper than that expected for a simple Keplerian rotation (IWeintroub et al.1 120081 ) . 



The radio emission of MWC349A is approxi mately prop ortional to u os , consistent with 



what is expected from an ionized stellar wind (jOlnonl 119751 ). High resolution VLA imag- 



(White & Becker 


1985; 


Cohen et al. 


1985). 


Tafova. Gomez, & RodrieuezI ( 


2004 


) argue that 


the dark lane with no free-free emission must be due to neutral gas. 


Tafova, Gomez, & Rodriguez 



show that the flux density increases with frequency as is' 



0.67±0.03 



and the angular size de- 



creases with frequency as v °- 74±0 - 03 ) confirming the presence of a biconical thermal wind that 
expands at constant velocity. The slow (~ 50 km s~ x ) bipolar wind is believed to originate 
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from photoevaporation of the neutral disk (jHamann &: Simon! Il986l; iLugo. Lizano fc Garay 



2004) . Analysis of VLA data over a 20 year time period (ITafoya. Gomez. &: Rodriguez 



2004 



Rodriguez. Gomez &: Tafoya 1 120071 ) suggest that the radio emission is fading by ~ 1.5% per 
decade. Th e free-free emission still dominate s at millimeter-wavelengths and perhaps even 
at 100 (Ltm (jHarvey. Thronson fc GatleyHl979l ). 



De ep Ha imaging of MWC 349 shows a faint shell to the north of the star with a radius 
of 75" (IMarston fc McColluml 120081 ) . This shell is seen to the east and north-east of the 
star in a MIPS 70 /im archive image, suggesting that MWC 349 A is illuminating a faint low 
surface brightness H II region. 

Our SCUBA observations at 850 /im and 450 /im show a point-like source with flux 
densities that are consistent with free-free emission (Table I2"al Figure ED- However, if we 
assume that the wind is launched at ~ 15 AU from the star (ITafoya. Gomez, fc Rodriguez 



20041 ). the wind should become optically thin at ~ 450 /im, and the free- free emission will then 
be proportional to u~ 0,1 . We may therefore start to pick up some dust emission shortward 
of 450 /im. We currently have no observations that enable us to measure the mass of dust 
and neutral gas in the disk. 



5.2. Embedded HAEBE Stars 



All embedded HAEBE stars are essentially "classical" HAEBEs, i.e., they illuminate 
reflection nebulae and are still located within the clouds in which they were formed. On 
average they are younger than the isolated HAEBEs and many drive bipolar outflows and 
can be deeply embedded. Their spectral classifications, and hence the ages, are also more 
uncertain, because some of them are very faint or invisible in the UV and visible parts of 
the spectrum. Some of them also have much higher accretion rates and veiling than isolated 
HAEBEs. Many of the embedded stars that we have observed are early to mid-B stars 
illuminating large reflection nebulae. With one or two exceptions they are non-detections, 
i.e., they no longer have observable disks in the sub-millimeter. These non-detections will 
be discussed in the next section (Section ED- 

Lkha 198 and V 376 Cas are late B and A or F stars at a distance of ~ 600 pc 
(Table [TJ. Both stars have been well studied at optical, near- infrared, far-infrared and sub- 
millimeter wavelengt hs. LkHa 198 has a deeply embedded inf rared companion, LkHa 198 B, 

100 



6" to the north (ILagage et al.l Il993l ; iKoresko et al.l Il997l ) with a luminosity of 



Lm, i.e., roughly similar to that of LkH a 198. LkHa 198 drives a bipolar molecular outflow 
(jCanto et al. 1984 ; Nakano et al. 1990 ) and excites several Herbig-Haro objects. A deeply 
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embedded, cold protostellar object, Lk Ha 198- mm, ~ 19" to the northwest, dominates the 
sub-millimeter emission in this region ( Sandell and Weintraub 1994al Henning et al. 1998h 
and may power the large bipolar molecular outflow ( ICanto et al.lll984t iNakano et al.lll990[ ) 
that originally was believed to originate from LkHa 198. LkHa 198 and LkHa 198 B lie 

inside a n elliptical limb-brightened nebulosit y. Bot h stars appear to ex cite Herbig-Haro 

objects (I Corcoran et al.l 1 19951 ; iKoresko et al.lll997l ). iKoresko et al.l (Il997| ). using adaptive 
optics compensated speckle observations in the H band, found a bar-like structure extending 
to 3" on each side of LkHa 198. They found LkHa 198 B to be extended at J and H as well , 
suggesting that both stars may be surrounded by circumstellar disks. iPerrin et al.l (120041 ) 
used laser guide star adaptive optics and near-infrared dual-channel imaging polarimetry to 
uncover a strongly polarized biconical nebula around LkHa 198 with a radius of 5", as well 
as a polarized, extremely blue, jet-like feature emanating from LkHa 198 B. Although these 
results strongly suggest that LkHa 198 has a disk, they still largely trace the surrounding 
infalling envelope. V376Cas has a bar feature, it illuminat es a sickle-shaped nebulosity, is 



very highly polarized and excites several Herbig-Haro knots ( iPiirola. Scaltriti fc Coyndll992 



Corcoran et al.lll995l ). All of this suggests that it is also surrounded by a circumstellar disk . 
Such a disk has been confirmed by speckle interferometry at H and K' (ISmith et al.ll2004j ). 
which show that V376Cas is surrounded by a compact, nearly edge-on flared disk with a 
radius of 100 mas. 

C O J = 2 — > 1 mapping of the central part of the CO outflow with 20" spatial reso- 
lution (jSandell and Weintraublll994al ) suggested that there is more than one outflow in the 
vicinity of LkHa 198 and V376Cas. LkHa 198 and/or LkHa 198 B could be driving one 
outflow, V376Cas another, and LkHa 198-mm could be the source for the large scale CO 
outflow. High spatial resolution CO J = 1 — > imaging (~ 6") with the Berkeley Illinois 
Mary land Array (BIMA) su pplemented by single dish mapping to fill in missing short spac- 
ings (IMatthews et al.l 120071 ) confirm that there are several outflows. V 376 Cas drives one 
outflow, LkHa 198-mm drives another, while there appear to be two outflows associated 
with LkHa 198. LkHa 198 definitely powers one outflow. The second outflow is closer to 
LkHa 198 B, but since it does not show any high velocity wings, they associate the driver of 
the outflow with a yet unseen source between LkHa 198 and LkHa 198 B. 

Our 850 /j,m map (Figure [7]) shows extended dust emission over ~ 3'. In the northeast, 
the emission becomes negative due to more extended and possibly stronger emission to 
the west, i.e., an artifact due to chopping. The 850 jttm dust emission is dominated by the 
protostellar source LkHa 198-mm. There is a tail of emission to the southeast from LkHa 198- 
mm and a fainter ridge extending north towards V376 Cas. The southwestern extension can 
be explained by an extended source, ~ 11'.'4 x 5'.'6, centered on LkHa 198 and its infrared 
companion with a position angle similar to that of the binary system (Table l2al) . The peak of 
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emission is close to LkHa 198 B, which agrees with the 1.3 mm map by lHenning et all (119981 ). 
which has slightly better angular resolution and shows a secondary peak on LkHa 198. Both 
the 1.3 mm map and the 850 /xm map suggest that LkHa 198 and LkHa 198B are surrounded 
by dust disks. We estimate that at most one third of the emission from the LkHa 198 binary 
is due to LkHa 198. If we assume a dust temperature of 50 K and a j3 = 1.5, we find that 
LkHa 198 has a disk mass of < 0.15 M Q . 



Interferometer observations ( IDi Francesco et all 1 199 71 ; iMatthews et alll20070 at 2.7 mm 



did not have enough sensitivity to detect LkHa 198 or V376 Cas, although lDi Francesco et al. 



did detect LkHa 198-mm. They found it to be unresolved or marginally resolved, while our 
SCUBA observations show it to be quite extended (FWHM ~ 20") with an integrated flux 
density of 1.78 Jy. We do see an enhancement in continuum level near V 376 Cas, but the 
position we deduce is off by 3" from the optical position and it could simply be due to emission 
from the surrounding dust cloud. However, observations at 450 or 350 /im at JCMT or CSO 
would have enough spatial resolution to resolve the HAEBE binary from the protostellar 
source and to determine how much of the emission originates from LkHa 198-IR and how 
much originates from LkHa 198. 

V 8 92 Tau is a heavily obscured (A y ~ 10 mag) late B star with a str o ng in frared 
excess (lEliaa 1 19781 ; IStrom & Stroml Il994l ). It was first discovered by lEliad (119781 ) in a 
near-infr ared survey of th e Taur us dark cloud complex. Inspection of H2 column density 
maps by iGoldsmith et all (120081 ) show it to be located in an opaque area of L 1495. Mul- 
tiepoch trigonomet ric parallax observations with the Very Long Baseline Array (VLBA) 
([Torres et all l2009f ) indicate that L 1495 is on the near side of the Taurus complex at a dis- 
tance of ~ 130 pc. iHernandez et all (12004 ) derived a spectral type of B8 Ve for V 892 Tau, 
while IStrom fc Stroml (I1994T) classified it as B9. Spectral types as late as A6 have been 
reported in the literature (IBerrilli et allll992l ). near-infrared speckle observatio ns and mid 



infra r ed interferometry show that it is a close bina ry with a separation of 50 mas (jSmith et al. 



20051 ; iMonnier et alll2008l ). iMonnier et all (120081 ) analyzed the SED for the binary and de- 
duced a bolometric luminosity of 350 L (corrected to 130 pc), with both stars having the 
same spectral type. The inferred luminosity is too low for two B8 V stars, an d the stars are 



more likely spectral type B9 V. Mid-infrared segment-tilting interferometry (IMonnier et al 



20081 ) shows a circumbinary disk inclined at 60° with an inner hole diameter of 33 AU. 



V892 Tau w as first detected in dust continuum at 1.3 mm by iBeckwith et all (Il990l ). 
Mannings! (119941 ) carried out extensive (sub-)millimeter photometry (1.3 and 1.1 mm, 850, 
800, 750, 450, and 350 /im) of V892 Tau with the single pixel detector UKT14 on JCMT, 
an d found it to have a v ery flat spectral index, a = 2.18 ± 0.10. It was mapped at 1 . 3 mm 
by iHenning et all (119981 ). who found it compact and unresolved. IDi Francesco et all (119971 ) 
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detected it in the dust continuum with the VLA at 7 mm and at 2.7 mm with the BIMA array, 
but di d not resolve the continuum emission. Recent CARMA observations by lHamidouche 
( 120101 ) at 2.7 and 1.3 mm resolve the disk with a size of 0'.'4 at 2.7 mm. 



V892 Tau is a very bright X-ray source (IStrom fc Stromlll994f). It is only seen in hard 



X-rays, suggesting th at it has low or no accretion (Hamid ouche et al.ll2008). How ever, it is a 



thermal radio source (ISkinner. Brown, fc Stewartlll993 



Di Francesco et al. 



19971 ). suggesting 



that the re is still some accreti on going on, although it is not known to drive a molecular 



outflow. IStrom fc Stroml (11994 ) deduced an age of ~ 3 My r by comparing its 



H-R diagram with two different evolutionary track models. iHernandez et al.l (120041) as sumed 



ocation on the 



it to be a B8 V star and found it to lie on the main sequence. IStrom &: Stroml (119941 ) found 
it somewhat puzzling that the star would appear older than nearby T Tauri stars, when 
it is still deeply emb edded in the molecular cloud core and illuminates a reflection nebula. 
Monnier et al.l (120081 ) find that it has an SED resembling a transition disk object and a clear 
inner gap in the disk, which is one of the main signatures of an evolved disk. However, since 
it is a binary, the gap in the disk is most likely caused by the binary. 

Our SCUBA images show that V892 Tau is unresolved at bot h 850 and 450 /mi (Fig- 
ure EJ), consistent with the small size found by iHamidouchd (|2010| ). We also retrieved and 
analyzed archive Spitzer MIPS images at 24 and 70 /im. At 24 /mi the star is heavily sat- 
urated and we could not deduce a flux density. The 70 /im image is of excellent quality 
and shows extended warm nebulosity east and north-west of the star (Figure [HD- We did psf 
photometry of this image and deduced a color corrected flux density at 70 /mi of 27.0 ± 0.4 
Jy. However, examination of the residuals from psf-fitting indicate that the star is mildly 
satura ted. The deduced flux is therefore strictly a lower limit, although probably not by 
much. |Pi Francesco et al.l (119981 ) mapped V 892 Tau with the Kuiper Airborne Observatory 
(KAO) at 50 and 100 /mi. They found much lower flux densities than IRAS. At 50 /mi they 
found a peak flux density of 59.1 ± 5.5 Jy and a source size of 25", resulting in an integrated 
flux of 86 ± 45 Jy. Although the star dominates at 50 /im, the KAO observations probably 
include some contribution from the extended warm nebulosity, which we see in the MIPS 70 
/im image. 

Isothermal graybody fits fit the millimeter and sub-millimeter data very well, but un- 
derestimate the flux at 100 /im and shorter wavelengths. Since the KAO beam at a 100 /tm 
is ~ 30", the observed 100 /mi flux density is an upper limit. The temperature of the dust 
emission is therefore somewhat uncertain; 40 K is probably a reasonable assumption. For an 
assumed dust temperature of 40 K we find a dust emissivity, j3 = 0.52, and a disk mass of 
0.009 M . If the dust temperature is 60 K we get = 0.40, and a dust mass of 0.004 M . 
In both cases the dust is somewhat optically thick at 850 /im, but less so for warmer dust. 
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As we can see from Figure [18] we can fit the millimeter and sub-millimeter data extremely 
well, but the 70 fim flux density is underestimated. The V892 Tau disk therefore appears 
very similar to the isolated HAEBE stars, which all require a warm dust component. The 
dust emissivity, which is reasonably well constrained, is rather low, suggesting that grain 
growth has already occurred in the disk. V 892 Tau has all the characteristics of an evolved 
disk. 

LkHa 101 is a heavily obscured (Ay ~ 10 mag), infrared bright Herbig Be star with a 
spectral type ~ BO Ve and lumi nosity of 4 x 10 4 L^, that is em bedded in the dark cloud 

L 1482 at a distance of ~ 700 pc (IHerbig. Andrews fc Dahirj|2004f) . It illuminates the reflec- 

tion nebula NGC 1579 and ionizes the H II region S 222 ( Herbig!ll960 . 1971 ; Herbig. Andrews &: Dahm 
20041 ) . It is the b rightest member of a deeply embedded c l uster of heavily reddened near- 



infrared sources ( Barsonv. Schombert and Kis-Halas] Il99ll; lAspin and Barsonyl [1994]) and 
faint cm-radio sources (feecker fe White! Il988l : Istine and O'Neall Il998h . most of which are 



probably young PMS sta rs. Interferometric imaging in the near- and mid-infrared (ITuthill. Monnier fc Dane 
20021 ; ITuthill et al.l 120021 ) shows that LkHa 101 surrounded by a dust disk with a hot inner 
wall facing the central star. The disk is seen nearly face on and shows a central hole or 
cavity. At 11.15 /im the size of the disk is ~ 0''06 (10 AU). They also find that the star is 
a binary system with a low luminosity companion ~ 0''18 to the northwest; this result has 



been independently confirmed by speckle interferometry (lAlvarez et al.ll2004l ). LkHa 101 is a 
strong ra dio source with a spectrum consistent with a spherically symmetric constant veloc- 
ity wind (IBecker fc White! Il988l ; iGibb &; Hoarell2007l and references therein). Although the 
star has been well studied in the centimeter regime, there is a lot of scatter in the free-free 
flux densities and size, even if we only look at VLA A-configuration data (Figure [12]) . This 
may be due to contamination from the surrounding H II region and because different array 
configurations resolve out more or less of the extended emission. It is also possible that the 
free- free emission is variable. If we use all high-spatial resolution observations (i.e., VLA 
A-array and MERLIN) available in the literatu re we get a spectral i ndex a = 0.71 ± 0.07, 
which is steeper than but consistent with what IGibb fc Hoard (120071 ) found, i.e., a = 0.61. 



In the sub-millimeter (Figure ITT]) . LkHa 101 is located at the southwestern boundary 
of an extended elliptical cloud with a s ize of 3' x V . Al t hough the map covers five of the 
deeply embedded VLA sources seen by lStine and O'Neall (I1998I ). only LkHa 101 has a sub- 
millimeter counterpart. At 850 /mi the star appears somewhat extended with a size of ~ 
6'.' 5 x 2'.'6 or ~ 2800 AU. This size is probably an over-estimate, since the mapping was 
done in scan-map mode with relatively poor sky transmission. The dust disk, however, is 
not very massive, since most of the emission at 850 /iin is due to free-free emission. In the 
low signal-to-noise 450 fxm map the emission is completely dominated by the surrounding 
dust cloud and we cannot derive a flux density for the star. After subtracting the free-free 
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emission we are left with ~ 180 mJy at 850 /mi. If we assume f3 = 1 and a dust temperature 
of 50 K we find that LkHct 101 could have a disk mass of ~ 0.005 M Q , which is similar to 
the disks of isolated Ae stars. 

AB Aur is a bright, optically visible Al Ve star located on the front side of the dark 



cloud L 1517 in the Taurus- A uriga dark cloud complex at a distance of 144 pc ( Ivan den Ancker et al. 
19981 ; iHernandez et al.ll2004j ) . It is ~ 3' away from the classic al T Tauri star SU Aur (spectral 



class G2 Hie), which is also surro unded by an accretion disk (IWeintraub. Sandell. fc Duncan 
19891 ; Andrews fc Williamsll2005l ). AB Aur was included in the lMeeus et al.l (120011 ) sample of 
isolated HAEBEs. However, AB Aur is not an isolated HAEBE star, since it is still located 



withi n the dark cloud in which it was formed. Its age, however, 4 ± 1 Myr ( IDeWarf et al. 



20031 ) is similar to some of the younge r isolated HAEBEs. N ear-infrared coronagraphic imag- 
ing with adaptive optics on Subaru (IFukagawa et al.l 120041 ) show extended emission as far 
as ~ 4" from the star and resolve the emission into a spiral-like structure with at least two 
spiral arms. 

AB Aur has been exte nsively studied in the (sub-) millimeter, both with single dish tele 



scopes and interferometers ( Mannings! 1994 ^ Mannings fe Sargent 1997 ; Pietu. Guilloteau fc Dutrey 



2005; Semenov et al 



free e mission ( IGiidel et al 



2005 



Corder et al 



1989 



2005 



Lin et al 



2006). The star is detected in free- 



Skinner. Brown, fc Stewartlll993l ; iRodriguez. Zapata. &; Ho 



20071 ) but is very faint, emitting only 0.2 ± 0.02 mJy at 3. 6 cm. High spatial resolution (0"6 



- 2") interferometr i c imaging at m illimeter wavelengths (jPietu. Guilloteau fc Dutrey! 12005 



Corder et al.l 120051 ; iLin et al.l 120061 ) in various CO transitions and their isotopomers show 
that disk is seen nearl y face-on with an inclination angle betw een 20° and 4 ° and appears 
to be non-Keplerian. iPietu. Guilloteau fc Dutreyl (120051 ) and ICorder et al.l ( 120051 ) resolve 
the continuum e mission and find it asymmetric, possibly tracing the spiral structure seen 
m near- infrared. ILin et al.l (120061 ). who observed AB Aur in 12 CO J=3 ->■ 2 with the SMA, 
confirm by tracing the inner spiral arm seen in the near-infrared in several velocity channels 
that the disk deviates from Keplerian rotation. Their 870 continuum observations show 
a central depress ion toward s the star, suggesting the presence of a large inner dust-free hole 
around the star. ILin et al.l suggest that both the central depression and the spiral structure 
with non-Keplerian moti ons could be caused by a giant planet forming in the disk. This 
view is also supported by lOppenheimer et al.l (120081 ) . who did adaptive optics coronagraphy 
and polarimetry of the AB Aur inner disk. They find an azimuthal clearing in the dust at a 
radius of ~ 102 AU along with a clearing of dust at ra dii inside this annulus, suggestin g the 
formation of at least one planetary body at ~ 100 AU. IRodriguez. Zapata, fc Hoi (120071 ) also 
detected faint radio emission from a source located about 40 AU from AB Aur, suggesting 
that the star could have a low mass stellar companion. 
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We f i nd AB Aur unresolved both at 850 /j,m and 450 fim (Figure [T3|) . even though 
Lin et al.l ( 120061 ) resolved it at 870 jum with a size of 3".l x l'/9. The star also appears 



compact in a MIPS 70 fim image, but the star is somewhat saturated and we could not 
retrieve a flux density. Our flux density measurements with SCUBA (Table I2"a|) agree within 
errors with our less accurate UKT 14 photometry (Table [2b]) . Figure [T41 shows a least-squares 
graybody fit to all reliable (sub-)millimeter and far-infrared data of AB Aur, resulting in a 
dust temperature of 66 K and a dust emissivity index, (3 = 1.05. If one assumes a normal 



gas-to-du st ratio, the mas s of t he disk is ~ . 006 



Mp). This disk mass is consistent with the 



results of iLin et al.l (I2006j) andlCorder et al.l (120051). but low er by a factor of three compared 
to the mass obtained by lPietu. Guilloteau &: Dutreyl ( 120051 ) . 



VY Mon is a late B to early A star at the southern edge of the bright reflection 
nebula IC446 at a distance of 800 to 900 pc (jCasey &: Harperl I199CH ) . The star is heavily 
reddened and variable wit h a bo lometric luminosity of ~ 870 L for a distance of 800 pc . 
The, de Winter and Perea ( 119941 ) assigned it a spectral class B8:e, while iMora et al.l (120011 ) 
classified it as A5 Vep. The observed luminosity appears too high for an A5 star, and we 
have therefore adopt the s pectral classificat ion by iThe. de Winter and Perezl . i.e., B8. The 
star is a faint radio source ( Felli et al"1ll982l ). VY Mon drives a high- velocity outflow, which 
we mapped in CO J = 2 — > 1 on Figure [15] shows high velocity emission from about 

—28 to +18 km s^ 1 . The outflow is very compact and appears to be seen pole on, since the 
blue a nd red-shifted outflow lobes are completely overlapping (Figure [16]). IWang &: Looney 



(I2007T ) analyzed near-infrared and Spitzer IRAC observations of IC446 and find that VY 
Mon is part of a compact group of 26 YSOs within a radius of 0.5 pc. 



VY Mon was mapped at 1.3 mm by IHenning et al.l (jl998l ). who found strong, compact 
1.3 mm emission centered on the star and faint extended emission from the dark cloud 
south of VY Mon. We have no SCUBA observations of VY Mon, and our single beam 
photometry with UKT14 only g oes down to 800 //m ( Table l2b]) . At 1.3 mm our observed 
flux density is higher than what IHenning et al.l (119981 ) measured from their higher spatial 
resolution IRAM observations, suggesting that we pick up some extended emission from 
the surrounding cloud in our larger beam. A graybod y fit to the (sub- ) millim eter data 
underestimates the far-infrared flux densities observed by ICasey fe Harperl ( jl990l ) and gives 
f3 — 1.9 and a disk mass of 1.9 M . This mass is almost certainly too high to be from the disk 
alone and suggests that our single dish observations may include a substantial fraction from 
the surrounding inf ailing envelope. The high value of j3 is also consistent with unevolved 
(small) dust grains, which is additional evidence that this flux is capturing emission from 



3 The original data for these CO observations are missing in the JCMT archive and in our personal archive, 
but we have preserved in our records and present in this paper the maps derived previously from those data. 
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the inf ailing envelope. 



AFGL 961 is a deeply embed ded (Ay > 18 mag), young, high-luminosity, double star 
system with a separation of ~ 6" (jCastelaz et al.lll9851 ). Both stars appear to be early B 



stars and ar e surrounded by extended nebulosity in the near-infrared ( jCastelaz et al.lll985 



Aspinl 119981 ) . AFGL 961 is southwest of the the H II region NGC2244 (the Rosette nebula) 
and part of the Mon OB2 association, which lies at a distance of 1.4 - 1.6 kpc ( jTurner 
19761 ; iHensberge. Pavloski fe Verschuerenl 120001 ) . Here we adopt 1.5 kpc for the distance to 
AFGL 9 61. The double star sys t em is located at the center of a large bipolar molecular 
outflow (ILada and Gautierlll982t iDent et al.l 120091 ). The main outflow has an extent of 5' 
at a P.A. of 45°, but IDent et al.l found at least one other, more compact outflow orthogonal 
to th e main outflow. Early millimeter continuum observations ( iHenning. Pfau fe Altenhofi 
19901 ) easily detected AFGL 961, but had insufficient spatial resolution to identify which star 
is responsible for t he dust emis s ion. T he same was true with a more recent (sub-) millimeter 
imaging survey by iKlein et al.l (120051 ). who showed that AFGL 961 is embedded in a dense 
sub-millimeter core with a size of 0.3 pc and a mass of ~ 550 M . 

The eastern component, AFGL 961 A (al so known as AFGL 961 E ) is more deeply em- 
bedded and a faint free-free emission source ( IBally fe Predmord 119831 ). with stronger ther- 
mal excess in the near-infrared t han AFGL 96 1 B (also known as AFGL 961 W), although 
both stars have Bra in emission (lAspinl I1998I ). AFGL 961 B shows the 2.3 ixm CO band- 
heads in emission, while the eastern source shows a featureless spectr um except for the 
Bra emission. Analysis o f 2MASS and deeper near-infrared imaging ( ILi fe Smith! 12005 



Roman-Zuniga et al.l 120081 ) shows that the majority of stars in the Rosette molecular cloud 
complex are contained in clusters. In the near- infrared AFGL 961 with its two early B stars, 
ap pears to be a very small clus ter that is barely detectable in the deep near-infrared imaging 
by lRoman-Zuniga et al.l ( 20081 ). who identify it as PL 06. Recent high-spati al resolution mid- 
infrared and sub-millimeter imaging with the SMA (jWilliams et al.l 120091 ) show that PL06 
is a deeply embedded, very young cluster. They find three continuum sources at 1.4 mm. 
The strongest one, SMA 1 coincides with AFGL 961 A, while two of them have no near- or 
mid-infrared counterparts. SMA 2 appears to be starless, but accreting, i.e., a pre-stellar, 
collapsing dense core, while SMA 3 drives an outflow, and is possibly driving the large NE- 
SW outflow. SMA 3 is therefore a Class protostar. They also identify a third mid-infrared 
source, AFGL 961 C, ~ 30" west of A, with strong mid-infrared excess, but they did not 
detect it in continuum at 1.4 mm. AFGL 961 A was also detected with the VLA at 7 mm, 



with a flux density of 8.1 ± 0.1 mJy (jZapata et al.l 120091 ). 



At 850 ixm we find a strong elliptical source centered slightly northwest of AFGL 961 A 
and superposed on relatively strong, extended emission from the surrounding cloud (Fig- 
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ure fT71) . At 450 /im the elliptical source is resolved into two compact sources. The eastern 
one is approximately centered on AFGL 961 A, possibly with some contribution from SMA 2. 
The second 450 /im source, AFGL 961 SMM overlaps with SMA 3, but is definitely centered 
east of SMA 3, suggesting that there are even more sub-millimeter sources in the AFGL 961 
cloud core. To learn more about this complicated young cluster, we therefore also retrieved 
and reduced Spitzer archive MIPS 24 and 70 /im images of AFGL 961. The 70 /im image is 
plotted in gray scale in Figure [TTl and overlaid with the 850 /im SCUBA image in contours. 
In Figure [T7| we also separately plot both the 70 and the 24 /im images over the same area 
mapped by SCUBA. AFGL 961 A and B are heavily saturated both at 24 and 70 /im, but 
SMA 3 coincides with a clear 70 /im peak. There is, however, a clear extension to the south- 
east towards A & B confirming that the western 450 /xm source has a 70 /im counterpart. 
What we see at 450 and 70 /im are probably two partially blended sources, one of which 
coincides with SMA 3. At 24 /im AFGL 961 C is also seen as a double source, but with the 
southea stern componen t stron ger than C. The second source (which we call 961 D) is also 



seen by I Williams et al.l (120091 ) at 20 /im as a faint extension to the southwest. At 70 /im 
there is no clear peak at C and D although both sources lie within the relatively bright 70 
/im core. A much stronger ridge of emission is seen protruding to the south-southeast of A 
& B. This ridge of emission is also seen in the SCUBA images, and suggests that there may 
be yet more deeply embedded sources hidden in the AFGL 961 cloud core. 

Analysis of the embedded cluster population is outside the scope of this paper, but 
AFGL 961 A has all the characteristics of a very young Herbig Be sta r, which most l i kely is 



surrounded by an accretion disk and an infalling envelope. Although IWilliams et al.l (120091 ) 
did not find an outflow from the star, the spectral index of the free-free emission between 
5 and 15 GHz is con sistent with the free- free emission originating in an ionized wind or jet 



(ISnell fc Ballylll986l ). A graybody fit to all (sub-)millimeter data gives a dust temperature 
of ~ 37 K, a dust emissivity, (3 ~ 1.51, resulting in a total mass of the disk/envelope of 12.5 
M . This fit overestimates slightly the flux density of SMAl at 1.4 mm, but is rather close 
to the sum of SMA 1 and 2, suggesting that we pick up emission from both sources in our 
much larger SCUBA beam. On the other hand the fit severely underestimates the 7 mm 
flux density, which only comes from A (SMAl), suggesting that part of the emission may 
be free-free emission, rather than thermal dust emission. 

To get a better idea of the true luminosity and mass of AFGL 961 A we have addi- 
tionally explored the archive of two-dimensional axisymmetric radiative transfer models of 
protostars calculated for a large range of protost e llar m asses, accretion rates, disk masses, 



and disk orientations created by iRobitaille et al.l (120061 ). This archive also provides an on 



line regression tool for selecting all the best fit SEDs. We use all near-, mid- infrared and 
millimeter photometry that resolve the star as well as a Spitzer IRS Short-High (SH, i.e., 9.9 



-28 - 



- 19.6 fim) spectrum from the Spitzer archive. Our results are similar to what I Williams et al. 
(120091 ) found using the same tool. The best fit models predict a stellar mass of 10 - 12 M , 
obscured by 12 - 18 mag of visual extinction. These models result in a disk mass of 0.02 - 
0.4 M , an envelope mass of 20 - 30 M , an envelope accretion ra te of 10~ 6 M m yr -1 , a 



stellar age of ~ 10 4 yr and a bolometric luminosity of ~ 3,000 L . IWilliams et al.l (120091 ) 
find a slightly higher luminosity, ~Lboi ~ 4000 L , but otherwise our results are very similar. 
De Wit et al.l (120091 ) modeled the observed SED of AFGL 961 A with a spherical model with 
a radial density profile described by a simple power law using partially the same data as 
we do. Their best fit model predicts a power law index of 0.5 and gives a total mass of 25 
M and a luminosity of ~ 6000 L . However, some of the data, like MIPS 70 /im and our 
SCUBA data included in their SED, includes contributions from the nearby sources SMA2 
and AFGL961B and therefore the derived luminosity is over-estimated. 

AFGL 961 A appears to be a very young star driving an ionized outflow (and possibly 
a mole cular outflow) t hat is surrounded by a large accretion disk. The modeling results by 
us and I Williams et al.l predict a protostellar mass of 10 - 12 M and a luminosity of 3000 - 
4000 L , which correspond to a B2 star. Since the star is still heavily accreting, it may end 
up as perhaps a Bl star by the time it reaches the main sequence. 

R Mon is a deeply embedded HAEBE star at the a pex of the NGC 2261 nebu l a, also 
known as Hubble's Variable Nebula. In the catalogue by iThe. de Winter and Perezl (119941 ) 
the spectral type is given as BOe. This is clearly incorrect , since the bolome tri c luminosity 
of the star is only ~ 740 L at a distance of 800 pc (jCohen et al.l Il985l ). iMora et al. 



(120011 ) classifie d it as B8 I 



determined by lNatta et al 



l ev, w hich is consistent with the stellar luminosity of ~ 450 L , 
( 19 931) . It is a b inary system with a low-mass companion located 



~ 0'.'7 to the northwest (jClose et al.lll997l ). The star drives a prominent molecular outflow 
and exci tes the HH object HH39, located in a small anonymous da rk cloud ~ 7' north of 
R Mon JCanto et all Il98ll : Ijones fc Herbkl Il982l : brueel et all Il984h. R Mon was the firs t 



HAEBE star to be detected with a millimeter-interferometer (jSargent fc Beckwithl 119871 ). 
They did not have enough sensitivity to detect the star with the OVRO three-antenna array 
in continuum at 3 mm, but detected 13 CO J = 1 — > in th e envelope surrounding the 



star. The star is a faint VLA source, 0.29 ± 0.03 mJy at 6 cm ( ISkinner. Brown. &: Stewart 
19931 ). with a spectral index in the radio regime of a = 0.9 ± 0.4, consistent with wind 



excited free-free emission. The firs t (sub-) mi ll imete r (1300, 1100, 800, 450, and 350 /im) 
detection of R Mon was reported by iManningd (119941 ). who als o determined t hat R Mon had 
a disk/envelope with a substantial mass of 0.24 M . Recently iFuente et al.l (120031 ) detected 
the continuum emission from R Mon at 2.7 mm using the PdB interferometer. They found 
a total flux of 6.4 mJy and claimed that the emission was extended on the 3" - 4" level. 
They also observed R Mon with the VLA at 1.3 and 0.7 cm, but do not report integrated 
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flux densities. In a follow-up paper (IFuente et al.ll2006l ). they observed R Mon with the PdB 
interferometer in the continuum at 2.7 and 1.3 mm and in the rotational lines of 12 CO J = 
1—^0 and J = 2 — > 1. Although they find a Keplerian signature in their CO images and 
derive a mass for the central star of 8 M , based on the bolometric luminosity R Mon is a 
late B (see above) and cannot be as massive as 8 M©. The derived CO rotation curve is 
most likely heavily contaminated by the strong outflow from the star. 

We have observed R Mon at various wavelengths between 350 and 1300 jttm in the dust 
continuum with the single pixel bolometer UKT14. The star was detected at all wavelengths 
(Table l2b|) . Peak-up obser yations of the star s uggest that the emission is compact, which 
agrees with the results by IFuente et al.l ( 120061 ). who barely resolved it at 1.3 mm with a 
size of 0'.'3. Fitting an isothermal graybody model provides a dust opacity index j3 = 1.3 
and a disk mas s ~0 .1M ff (Table H]). The disk mass we derive is much larger than found by 



Fuente et al.l (120061) . because we find a much higher dust emissivity, 1.3 vs 0.3 - 0.5 deduced 



by IFuente et al.l Our observations cover a much larger frequency range than IFuente et al. 
and, since the free-free emission is negligible in the sub-millimeter regime, our deduced dust 
emissivity is rather well constrained (see Figure [T^]) . It is possible, however, that the envelope 
may contribute to the observed dust emissivity, especially since most of our data come from 
single pixel observations. 

HR5999 (HD 144568) illuminates a bright refle ction nebula in the c entral part of the 
Lupus 3 dark cloud complex at a distance of 208 p c (IPreibisch et al.l 120071 ). The star varies 
strongly and irregularly on time scales of weeks ( Perez et al. 1992h . H R5999 was classi- 



20071 ). Ivan Boekel et al.l (120051 ) found a similar spectral type as lBlondel fc Tjin A Djiel . A5 



fled as A6 III with an error of about one subclass ( iBlondel fc Tjin A Djiel 120061) and forms 
a ~ 4 5" wide proper-motion binary system with the Al.5 st ar HR6000 (Preibi sch et al. 



7 III/IVe, and deduced a stel 
a strong far-infra red excess 



ar mass of 3.2 ± 0.5 M ^ and an age of 0.5 Myr. The star has 



Hillenbrand et al.lll992l ) and shows broad silicate emission in 



the mid-infrared ( Siebenmorgen et al. 200o|). but it is a rather faint millim eter source and it 



was not detected at 1.3 mm ( IHenning et al.lll994j ). IPreibisch et al.l ( 120071 ) observed the star 
in the mid-infrared with MIDI at the ESO Very Large Telescope Interferometer (VLTI) and 
was able to resolve the circumstellar disk around the star. They found that the disk is highly 
inclined, at 58°, and truncated at 2.6 AU, likely due to a close binary. We barely detect 
the disk as an unresolved point sou rce at 850 jttm (Table l2a|) . A graybody fit to the 1.3 mm 
upper limit by lHenning et al.l ( 11994] ). our 850 fim flux density, the 100 fxm IRAS upper limit, 
and the IRAS 60 /im flux density, yields a dust temperature of 41 K, a (3 = 1.75, and a dust 
mass of ~ 0.009 M , which is comparable to the mass estimated by ISiebenmorgen et al. 



( l2000l ) 
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MWC 297 is a highly reddened (Ay ~ 8 mag) early t ype star of spectr al class B1.5 
Ve, with a luminosity of 3 x 10 3 L Q and a distance of 250 pc (IDrew et al.lll997l) . It has very 
strong and variable Balmer line emission 



Finkenze 



strong radio source (jSkinner. Brown. & Stewart 



1993 



ler fc Mundt 



Drew et al 



1984) and is a relatively 



19971 ). |Pi Francesco et al. 



(119981 ) found the far-infrared emission from MWC 297 to be extended with a size of about 
one arcminute at 100 fim. However, at 1.3 mm the map of MWC 297 ( jHenning et al.lll998[ ) 
shows that the millimeter emission from the star appears to be compact, and that there is 
strong emission from the surrounding dust cloud to the northeast of the cloud. The 1.3 mm 
map also shows compact emission from a very red 2MASS s ource (J18273709-03 4938), with 
a K magnitude of 8, located about 40" west of MWC 297. iManoi et all (boosl ). using the 
SMA, found the 1.3 mm continuum emission to be unresolved with a size of < 80 AU. They 
did not detect molecular gas emission in the J=2 — > 1 rotational transitions of 12 CO, 13 CO, 
and C 18 0. They deduced a very low dust emissivity index, (3 = 0.1 - 0.3, and suggested that 



the lo w f3 index could be due to grain growth if the dust is optically thin. lAlonso-Albi et al. 



(120091 ) observed the star with the PdB interferometer at 1.3 and 2.6 mm and with the VLA 
at 1.3 cm and 7 mm. They also found the millimeter emission to be compact and modeled 
the SED with a disk with a radius of ~ 28 AU seen approxim ately face-on, consistent with 



20081 ). who found that the disk is 



near-infrared interferometric observations (lAcke et al. 
seen with an inclination of < 40°. 

Our sub-millimeter maps (Figure 120]) agree with the 1.3 mm map. MWC 297 dominates 
the emission at 850 /xm, while the star is quite faint at 450 /im. At 450 fim the emission 
from the surrounding cloud is much stronger than from MWC 297, suggesting that the sub- 
millimeter SED is exceptionally flat. Follow-up observations by Sandell et al. (in prep) with 
BIMA and the VLA, show that the observed shallow SED at millimeter wavelengths is due to 
free-free emission. They observed MWC 297 with th e VLA in the AnB configura t ion at 6 and 
2.0 cm and also re-analyzed VLA data obtained by lSkinner. Brown, fc Stewartl (119931 ) . The 
Skinner. Brown, fc Stewartl CnD data show that MWC 297 powers a large bipolar ionized 
jet. A fit to the VLA data shows that the radio emission has a spectral index, a ~ 0.97, 
while the size decreases as z/ -0 8 , as is typically seen in collimated thermal jets. If we include 
data at wavelengths out to 3 mm, we get a spectral index, a ~ 1.03, i.e., essentially the 
same as we got from the more uncertain fit from VLA data alone (Figure l2Tj) . The thermal 
emission from the jet can therefore explain all the emission at 3 mm and there is no evidence 
for any excess due to cold dust emission. There is a hint of a dust excess at 850 fim and 
almost certainly some excess at 450 /iin, suggesting that the wind becomes optically thin 
somewhere around 300 - 600 GHz. MWC 297 looks therefore very similar to MWC 349, 
in that the emission in the millimeter and sub-millimeter regime is dominated by free-free 
emission. The bipolar morphology seen in the VLA maps suggests that the disk is seen 
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nearly edge-on, as is the case for MWC349. Both stars appear to be surrounded by disks, 
which are largely ionized and only contain a small fraction of neutral gas and dust. 

The BD+40°4124 r egion f o rms a small pre-main-sequence stellar cluster with three 
bright emission line stars. iHerbigl (119601 ) identified the first two, BD+40°4124 and V1686 
Cyg, as HAEBE stars. The third star, V1318 Cyg, is a dee ply embedded binary system 
with a separation of ~ 5" and a total luminosity of ~ 1600 L ( lAspin. Sandell fc Weintraub 
1994J ). The southern member of this bi nary, V1318 CygS, is more d eeply embedded with 
a visual extinction of 10 mag or more. lAspin. Sandell fc Weintraub! ( 119941 ) found that the 
position of V 1318 CygS coincides with the position of a compact sub-millimeter source that 
has a luminosity similar to that of a HAEBE star. None of th e stars in the BD+40°4124 



cluster show any free- free emission at centimeter wavelengths CSki nner. Brown, fc Stewart 



1993l : IPalla et alll995f ). except perhaps BD+40°4124, near which 
detected faint emission at 3.6 cm and 6 cm. 



Skinner. Brown, fc Stewart 



Palla et al.l (119951 ) discovered a com pact molecula r outfl ow and an H2O maser, which 
appear to be excited by V1318 CygS. lLooney et al.l (120061 ) . however, who observed the 



region in the 3.1 mm continuum with BIMA, found that the dust emission is offset by l'.'l 
to the north-east from V1318 CygS, and argue that the dust emission originates from a 
deeply embedded intermediate mass pro tostar V1318 CygS-mm, rather th an from either 
component of V 1318 CygS identified bv lAspin. Sandell fc Weintraub! (119941) . These results 
are supported by the less accurate results by Di Francesco et al. ( 19971 ) . who detected the 
star in continuum at 2.7 mm with the PdB interferometer, and found a similar offset, 0'.'9, 
from the southern co mp onent . M ore importantly, VLB A observations of the H2O masers 
toward BD+40°4124 (jMarvell 120051 ) indicate that it is associated with the sub-millimeter 
source rather than with the optically visible star. It is therefore the invisible millimeter 
source that drives the outflow, not V1318 CygS. 



Figure E2 shows the 850 and 450 /im SCUBA images of the BD+40°4124 region. It 
is clear from these images that neither BD+40°4124 and V 1686 Cyg are associated with 
any significant sub-millimeter continuum emission. Because of the extended emission from 
the molecular cloud in which these stars were born, we cannot place a very stringent limit 
on the dust emission; at 850 /im our upper limit is ~ 35 mJy beam -1 . Examination of 
MIPS 24 and 70 /im images show that all HAEBE stars in this region are saturated at 24 
/im (Table [3]). At 70 /im there is some extended emission near BD+40°4124, most likely 
heated dust from the surrounding reflection nebulosity. There is strong, compact, 70 /im 
emission from the vicinity of V 1318 Cyg S (Table [3]), almost certainly coming from the deeply 
embedded protostar, V1318 CygS-mm. This is where we find a strong, extended (8" x 3") 
sub-millimeter source (Table l2a|) . Although it is possible that some of the emission could 
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be due to V1318Cyg S, we assume that all the emission comes from the deeply embedded 
millimeter source. A graybody fit to the SCUBA and MIPS data give a dust emissivity of 
0.85 and a total mass of 1.9 M Q (Table Hj), suggesting that we see the disk and envelope of 
a young high-mass protostar. 

PV Cep is an irregularly variable star at the head of the cometary n ebula GM29 



that illuminates this fan-shaped nebula; it is located at a distance of 500 pc (ICohen et al. 



981 ). Its spectral class ifi cation is uncertain, because the optical spectrum is strongly veiled 



(Hernandez et al. 


2004 


)• 


Cohen et al. 


(1981) 



perhaps A5. T he latter would be consistent with t 
~ 80 - 175 L fevans II. Levreault fc Harvev!ll986 : 



re observed bolomet ric luminosity of 



Abraham et al.l 120001) . It drives a gi- 



ant p arsec-scale Herbig-Haro outflow (IReipurth. Bally, fc Devinelll997l ) and molecular out- 
flow lLevreaultlll984l: lArce and Goodman! l2002al Jbl). The outflow appears to be precessing 
( Arce and Goodman! 2002bh . sugg esting that PV Cep co uld be a binary system, although no 
companion has yet been detected. iTorrelles et al.l (119861 ) discovered an H2O maser 5" to the 
south of PV Cep, which could mark the location of another young star, because H2O masers 
excited by low-luminosity stars are usu ally close to t he star. More recent observations show 
H 2 maser activi ty close to PV Cep (IMarvell 120051 ). PV Cep coincides with a faint VLA 
source at 3.5 cm (lAnglada et al.lll992l ). 



PV Cep was mapped at 1.3 mm by (IFuente et al.lll998af). who found a c ompact source 
centered on the star. It was recently image on CARMA by iHamidouchd (120101 ) . who resolved 
it at 2.7 and 1.3 mm with a size of ~ 0'.'5 (250 AU). Our ma ps at 850 11m and 45 /im 
(Figure [23]) go deeper than the single dish map at 1.3 mm by IFuente et al.l (Il998al ). We 
resolve PV Cep in the sub-millimeter with a size of 1200 AU x 700 AU at PA 74° (Table EaJ, 
which is roughly perpendicular to to the outflow. We also see faint, extended emission with 
a diameter of ~ 1' from the dense surroundin g cloud core at bo th 850 and 450 fim. PV Cep 
was found to be unresolved at 50 and 100 /im (INatta et al.lll993l ). which is in good agreement 
with the compact size that we find in the sub- m illimeter. Our observed source size, however, 
is much larger t han what IHamidouchd (120101 ) found with CARMA and the 1.3 mm flux 
density found by IFuente et al.l (Il998al ) is about twice as large as what is seen by CARMA, 
suggesting that the array filters out so me of the extended em ission. A graybody fit to the 
(sub-) millimeter single dish data, ISO ( Abraham et ali 2000 ). and KAO far-infrared data 
(Figure [24]) gives a dust temperature of 29 K, a dust emissivity index, (3 = 1.59, resulting 
in a mass of ~ 1.0 M . This mass estimate is much larger than what one would expect to 
find even for a young circumstellar disk. It therefore see ms that the s i ngle d ish data maybe 
dominated by the compact infalling envelope, although IHamidouchd (120101 ) finds 0.8 M , 
suggesting that even array observations are dominated by the envelope. 
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V 645 Cyg powers a low velocity molecular outflow in the north-south direction ( jSchulz et al. 



1983 ; IVerdes- Montenegro et al 



199lh 



The star coincides with a f aint, somewhat exte nded 
VLA source ( Girart et al. 2002 ) roughly aligned with the outflow. Girart et al. (j2002 ) also 
found a faint extended, ~ 6' x 2', H II region close to the star, which c ould be excited by 
V645 Cyg. Since V645 Cyg is know n to power a 1665 MHz OH maser (IMorris and Kazes 
19821 ). and a Class II CH 3 OH maser (jSlysh et al.lll999l). it must b e a hig h- mass star, since 
these masers are only excited by high mass stars. IClarke et al.l (120061 ) classified it as a 
transition object between a massive YSO and a young Oe type star in a weak H II region. 



Di Francesco et al.l (119971 ) marginally detected V645 Cyg at 2.7 mm with the IRAM 
array. We find V 645 Cyg to coincide with a bright sub-millimeter source embedded in narrow 
east-west cloud ridge which turns mostly northward east of the star. The emission is clearly 
extended (Table I2"a|) . but at a distance of 4.2 kpc it is clear that the sub-millimeter emission 
is dominated by the dense envelope, rather than a circumstellar disk. An isothermal fit to 
the data in Table l2al gives a mass of 85 M for the elliptical sub-millimeter source, which is 
far more than what one would expect even for a very young circumstellar disk. 

MWC 1080 A must be an early B star, because it has strong B almer and Fe II lines 
showing P Cygni type emission pr ofiles (Finkenzeller fe Mundtl 1 19841 ) . It drives a promi- 



nent bipolar outflow seen pole on (ICanto et al 



984). Distance estimates vary. Here we 



have adopted 2.2 kpc (Tabled]) from ICanto et al.l (ll9841 ). MWC 1080 is a hierarchical mul- 
tiple system. The prim ary, MWC 1080 A, is an eclipsing binary. The period of 2.9 days 



(jShevchenko et al.lll994f). MWC 1080 B is also a PMS star and mo re heavily reddened than 



1997 



Po 



the primary (ILeinert. Ri chichi fe Haas 
found ~ 5"2 east of the primary (jPolomski et al 



omski et al.ll2002f). A th ird companion is 
120021: IPirzkal et al.lll997f). MWC 1080 i 



strong far-infrared source with a source size of ~ 5 0" - 60" (jEvans II. Levreault fe Harvey 



is a 



19861 : bi Francesco et alll998l : lAbraham et~al1l200ol ). VLA observations JGirart et al.lbooi ) 



show three free-free emission sources in the vicinity of MWC 1080 with VLA 4 coinciding 



Skinner. Brown. & Stewart 


(1998a 


) and 


Hennine et ali 



1993) at 3.6 cm with a flux density of 0.18 mJy iFuente et al. 



and found strong extended continuum emission towards MWC 1080 A and B. 

Our short integration 850 /im map (Figure 126]) is dominated by strong extended emission 
from the reflection nebula surrounding MWC 1080 and perhaps nearby embedded sources. 
MWC 1080 A, i.e., VLA 4, coincides with a faint emission peak in an extended east-west ridge. 
The other two VLA sources lack sub-millimeter counterparts. We consider MWC 1080 A a 
marginal detection with a flux density of ~ 250 ± 50 mJy. Thi s appears plausible since 
deep IRAM interferometer observations at 2.7 mm and 1.4 mm (jFuente et al.l 12003] ) show 
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extended emission towards MWC1080A. At 1.4 mm iFuente et all find that MWC1080A 
coincides wit h a fain t unresolved source. The observed flux density is much less than what 



Fuente et al.l (ll998al ) fouund from their single dish observations, suggesting that interferom- 



eter has resolved out most of the emission. If we assume that the observed 850 /im flux 
density comes from a circumstellar disk surrounding MWC 1080 A we can make some simple 
assumptions to estimate the mass of the disk. For a dust temperature of 50 K and a typical 
dust emissivity, j3 = 1.5, we get a disk mass of 0.009 M . 



6. Non-detections 



Most of our non-detections are early to mid-B stars. These stars illuminate large re- 
flection nebulae and are strong 60 and 100 /im sources in the IRAS point source catalogue. 
Only a few, however, appear to be surrounded by dusty disks that have enough cold dust to 
be detected with a single dish telescope in the sub-millimeter. Our sample also includes two 
stars which are not HAEBE stars but instead appear to be supergiants. Neither of these 
two, HD 316285 and MWC 300, are associated with any dust emission. These are discussed 
at the end of this section. 



MWC 137 is a BOep star that excites the H II region S266 flFichl Il993f ). Both the 
radio and the optical appearance are very similar: a ring-like nebulosity with a radius of ~ 
30" and two spiral-like arcs connecting to the centr al star (Figure |27|). A lthough MWC 137 
has long been known to be a radio source, see e.g., lAltenhoff et al.l (119761 ) . early single dish 
observations picked up emission from the exte nded H II region, not the central star. Radio 
emission from MWC 137 was first detected by lSkinner. Brown, fc Stewart! (119931 ) at 3.6 cm. 
They found the emission to be compact and unresolved. 

The star appeared to have b een detected at 1.3 mm in photometric observations at the 
CS O bv dHillenbrand et alll992h. It was mapped w ith the IRAM 30 m telescope at 1.3 mm 
by IFuente et al.l (Il998al ) and iHenning et al.l (119981 ) . Both maps show an extended narrow 
dust ridge extending over the nebulosity with a stronger emission peak ~ 30" northeast 
of the star and a fainter one to the southwest, both peaks approximately coinciding with 
the boundary of the H II region. There is only faint emission towards the star itself. Our 
850 /im map looks very similar, although we see a second arc bounding the nebula to the 
south. We see no enhancement in dust emission towa rd MWC 137; if any t hing the emission 
is close to a minimum at the position of the star. |Pi Francesco et al.l (119971 ) obtained a 
marginal detect ion of the star at 2 .7 mm with the IRAM interferometer, while the deep 
observat ions bv IFuente et al.l (120031 ) yield clear detections at both 2.7 mm and 1.4 mm. At 
1.4 mm IFuente et al.l derive a source size of 1"8 x 0"8 and find a total flux density of 10 
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mJy. iFuente et al.l (120031 ) also detected the star with the VLA at 1.3 cm and 7 mm, where 
they find it to be compact and unresolved. The VLA data suggest that MWC 137 has a 
thermal wind with a spectral index a = 0.8 - 1.0, i.e. similar to what we find for MWC 297 
and LkHa 101. It therefore appears that the free-free emission completely dominates at 
millimeter wavelengths. MWC 137 appears similar to MWC 349, MWC 297, and LkHa 101, 
but the free-free emission is much fainter, suggesting that it is more evolved and has already 
photo-ionized most of its disk. 



LkHa 215 illuminates the reflection nebula NGC 2245. It is a late B star. IValenti. Fallon fc Johns-Kru 



(2003) classify it as a B8 star based on low- resolution IUE spectra, while IHernandez et al. 



(j2004j) give it a spectral type B6 with an uncertainty of 2.5 s ubclasses. It is p art of 
the yo ung Mon Rl asso ciation, which has a distance of 800 pc (IHerbst et al.lll982f ). Al- 



though Mannings] (1994 ) report ed a detection of LkHa 215 with sub-millimeter photometry 
on JCMT, Fuente et al.l ( 1998a ). who mapped LkHa 215 and the surrounding reflection neb- 
ula in 13 CO J = 1 — > 0, CS J = 3 — > 2 and 1.3 mm continuum with the IRAM 30 m telescope, 
did not detect it. They found no CS emission in the vicinity of LkHa 215 and even though 
they detected 1.3 mm dust emission, it did not peak on the star. 

Figure [281 shows our 850 /im image overlaid on a POSS image of the region. One can see a 
faint dust ridge near LkHa 215. This emission, which is stronger east of the star, is probably 
not associated with LkHa 215, although we cannot exclude the possibility that LkHa 215 
is still surrounded by a small accretion disk. However, we find a much stronger unre- 
solved source, LkHa215-SMM, northwest of the reflection nebula (a2ooo.o = 06 ft 32 m 355265, 
^2000.0 — +10° 10' 18'/3), which has no optical or near-infrared counterpart. The flux density 
of this source is 152 ± 22 mJy at 850 /mi, 220 ± 80 mJy at 750 /im, and 875 ± 250 mJy at 
450 /im. 

HD 259431 is a B6 star that illuminates the reflection nebula NGC 2247. This star 
appears similar to LkHa 215. They are both part of the Mon Rl association and probably 
of similar age, 0.05 - 0.1 Myr (jTesti. Palla fc Nattalfl99~8l ). 



Mannings! (119941 ) reported a marginal 450 /im detection, but as is evident from the 850 
/im image (Figure [2H]) there is extended dust emission towards NGC 2247, but no emission 



associated with the position of the star. There are two faint 850 /im continuum sources 
outside the nebula; NGC 2247-SMM1: a 20 oo.o = 06^ 33 m 07599, <5 200 o.o = +10° 18' 19", and 
NGC2247-SMM2: a 200 o.o = 06 h 27 m 07576, <5 200 o.o = +10° 20' 22", with integrated flux 
densities of 0.38 mJy and 0.16 mJy, respectively. Neither source has an optical or near- 
infrared counterpart. 



LkHa 25, spectral type B7 (IHernandez et al.l 120041 ) . does not illuminate a reflection 
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nebulosity ( IHerbigl I196CH ) . nor was it detected by IRAS. Spitzer IRAC archive images show 
that the star is inside a large, roughly spherical nebula with a diameter of 6' - 8'. Although 
the star is not centered in the nebula, it is ~ 2' from the southern part of the nebula (which 
is also the brightest part of the nebula at 8 jum), suggesting that it may be t he illuminating 



source of the nebula. LkH a 25 has been observed by two groups at 1.3 mm (IHenning et al. 



19981 ; iFuente et al.lll998al ). but neither of them detected the star in millimeter continuum 



emission. We do not detect it either. Our sub-millimeter map (not shown) shows that the 
star lies in a region free of dust emission, although there is a small dust emission region less 
than one arcmin to the southeast of the star and a large dust cloud (~ 10' x 3') about two 
arcminutes west of the star extending mostly to the north. This cloud complex is almost 
certainly associated with NGC 2264. 

CD— 42° 1172 is a well-studied, bright BOep star illuminating a clump y, arc-shaped 

nebula with a size of ~ 40", showing both reflection and emission characteristics (IHutsemekers fc Van Prom 



19901 ). Although several studies have questioned whether CD— 42° 1172 is a PMS star (see 
e.g., BorgesFernandes et al. 2007), there should not be any doubt that this is a young star. 
It is embed ded in an opaque dark cloud and is the brightest member of a small cluster of 



PMS stars Jrlabart etaDl2003l: 



Boersma et al 



20091 ). Distance estimates vary. Here we have 



adopted 1.15 kpc ( Borges Fernandes et al. 20071 ). which is based on the interstellar absorp- 



tion lines of Na I and Ca II K lines. The star is associated with st rong far-infrared emission 
( Weintraub 1990 ; Natta et al. 1993 ) and sub-millimeter emission ( Mannings 1 1994), but it is 
questionable whether the emission comes from a circumstellar disk. iBoersma et al.l (120091 ) 
present an IRAC color image of the CD— 42° 1172 region. The star is saturated in all IRAC 
images and even a secondary peak north of the star is an artifact due to the heavily satu- 
rated CD— 42° 1172. The IRAC image shows that the star is surrounded by a one arcminute 
nebulosity glowing in PAH emission and it also shows a smaller nebulous region about 45" 
to the northwest. 



Henning et al.l (119981 ) mapped CD— 42° 1172 in 1.3 mm continuum with the IRAM 



30 m telescope. They found the star to be associated with an extended dust cloud, but 
there is no enhancement in the dust emission at the position of the star. Our 850 and 450 
/im images (Figure I3"0l look very similar to the 1.3 mm map and show no enhancement 
of the dust emission towards the star, only emission from the cloud ridge in which it is 
embedded. We see two sub-millimeter sources to the northwest of CoD— 42° 1172, neither of 



which have any near- or mid- infrared counterparts (IHabart et al.l 120031 ) . SMM1 is a point- 
like source offset — 6"8,+7"3 from the star in RA and Dec, i.e., a 2 ooo.o = 16 h 59 m 065286, 
£2000.0 = -42° 42' 00"7, with a flux density of 670 mJy at 850 /im. SMM2 is a strong (2.4 
Jy), extended source about 50" northwest of the star: 0:2000.0 = 16 fc 59 m 035911, £200 0.0 = 
—42° 41' 32"6. The latter is clearly seen in the 1.3 mm map by iHenning et al.l (119981 ) and 
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appears to coincide w ith the northwestern nebulosity seen in the IRAC image shown by 
Boersma et aD fl2009[ ). 



We therefore also downloaded and analyzed MIPS images of CD— 42° 1172, which has 
been observed as part of the MIPS GAL legacy project. The 24 and 70 /mi MIPS im- 
ages are shown in Figure [30j on which we we overlaid the 850 /tm SCUBA image with 
contours. In the 24 /zm image the nebulosity surrounding CD— 42° 1172 is completely sat- 
urated. One plausible explanation to why the nebulosity would be saturated is that it is 
an H II region ionized by the star or by a t hermal wind from the star. The star, how- 
ever, was not detected in the VLA survey by ISkinner. Brown, fc Stewart! ( 119931 ). nor was 
it seen by iThompson. Urguhart &: White! (120041 ). who used the Australia Compact Array. 
Thompson. Urguhart fc White! detected two sources at 13 cm. One is probably extragalactic; 
the other one, SF0 85a, is close to, but not coincident with the extended SMM2 northwest 
of CD— 42° 1172. The emission is also extended in the MIPS 70 /xm image, suggesting that 
this region may contain several embedded stars. 

TY CrA and HD 176386 illuminate the reflection nebulae NGC6727 and NGC6726 
in the Corona Australis m olecular cloud. TY CrA is a triple system with a close eclipsing 



binary (jCasey et al. 



an age of ~ 3 Myr (jCasey et al. 



1995) . The prima ry is a B8-B9V star with a luminosity of 67 L Q , and 



19981 ). HD 176386 i s a binary star wit h the primary being 



an A0V star and the secondary a low-mass K7 star (IMeyer et al.l 120091 ) 



Figure ED sho ws an 850 /im sub-image of a larg e SCU BA scan map, which has also 
been published by iNutter. Ward-Thompson fc Andre (2005) . It is clear from this 850 jj,m 
map that neither star has any sub-millimeter excess. iHillenbrand et al.l (119921 ) reported a 
1.3 mm detection of TY CrA, but as we can see from Figure |31j the emission observed by 
Hillenbrand et al. is from the dust cloud in which the star is embedded. Neither TY CrA 



nor HD 176386 show any near-infrared excess (IBibo et al.l 119921 ) . suggesting that they have 
completely dispersed their accretion disks. The nearby Herbig Ae/Be stars, R CrA and T 
CrA, about 5' to the southwest, appear to be much younger, and even though ne ither of 
them have been detected in the sub- millimeter (see Figure 2 in iGroppi et al.l ( 120071 )). they 
do show clear excess emission in the near- and m id-infrared, suggesting that they are still 
surrounded with accretion disks ( IBibo et al.lll992l ). 



HD 200775 is a B3 star illuminating the reflection nebula NGC 7023. We do not see any 
dust emission towards HD 200775 in NGC 7023 to a limit of ~ 25 mJy beam -1 (Table [2c]) . 
NGC 7023 was observed both in scan map and jiggle map modes. There is strong dust 
emission from the PDR region just north and northwest of the star and fainter emission to 
the east and southwest at the bounda r y betw een the reflection nebula and the surrounding 
dust cloud (Figure |32|) . iFuente et al.l (Il998bl ) show that the cavity around HD 200775 has 
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very likely been excavated by an energetic, bipolar outflow in an earlier evolutionary stage of 
the star, but that the star no longer shows any evidence of a bipolar outflow. They detected 
faint continuum emission at 3.4 mm with the IRAM interferometer, but the observed flux 
density, 3 mJy is abo ut what one would expect from a thermal wind. Recently, however, 
Okamoto etHI fj2009h detected HD 200775 with the SMA at 350 GHz with a flux density of 
35 mJy, i.e., consistent with our 3-a upper limit, see Figure [321 They also found unresolved 
excess emission in the mid-infrared using the 8.2 m Subaru telescope, confirming that the 
star is still surrounded by a disk. 



HD 316285 is listed as an extreme emission line object in IThe. de Winter and Perez 



( 119941 ); however, HD 316285 is a P Cygni t ype supergiant, po ssibly a luminous blue variable, 
and definitely not a young Herbig Be star (jHillier et al.lll998l ). It is extremely luminous and 
lies at a distance of ~ 1.85 kpc (Tabl e ED). It has very little or no hot du st, but shows strong 
far-infrared emission from cold dust jMcGregor. Hvland fc Hillierlll988h . 



van der Veen et al.l (119941 ) did sub-millimeter photometry with the single pixel bolome- 



ter UKT14 on JCMT and reported flux densities at 1.1 mm, 800 and 450 /im although the 
flux density at 800 /im was suspiciously low. In Figure [331 we present large 850 and 450 /zm 
images in the direction of HD 316285. Although there is strong sub-millimeter emission in 
the field, HD 316285, is in an area with relatively weak dust emission. There is a marginal 
peak in the dust emission at 850 /im, but not at 450 /im (Figure [33]) . We count it as a 
non-detection. The upper limit we derive from these observations is rather high, 70 mJy at 
850 /im for an unresolved source (Table l2"cl). 



MWC 300 was listed as a HAE BE star by iHerbigi (Jl960j) and included in the catalogue 



by IThe. de Winter and Perez 



or supergiant (lAUen fc Swings 



1994). However, it ha s also been classified as a B|e ] star 



19761 ; lAppenzellerl 119771 ) . iMiroshnichenko et al.l (120041 ) , who 



analyzed more than ten years of optical and near-infrared spectroscopy and mid-infrared 
imaging, were able to determine the distance (1.8 ± 0.2 kpc) and hence the luminosity (L = 
1.3 x 10 5 L Q ), which together with all the other observed characteristics and with detailed 
modeling identify this star as a B[e] supergiant and most likely a binary. 



Henning et al.l (119941 ) reported an approximately three-sigma detection of the star in 



continuum emission at 1.3 mm. We find no emission at all towards MWC 300 with a rather 



stringent upper limit, 7 mJy beam 1 (Table Be]) , suggesting that the results by lHenning et al. 
were spurious. 
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Discussion 



In this paper we present results from an extensive sub-millimeter survey of HAEBE 
stars, many of which have been detected both at 850 and 450 /im. First we note that our 
study includes a total of 39 objects that were potential HAEBE stars. However, HD 316285 
and MWC 300, which were classified as potential HAEBE stars when we surveyed them, are 
now classified as supergiants. We therefore have observed a total of 37 objects now thought 
to be HAEBE stars. Of these we have detected 27 stars (73%) in emission at 850 fj,m or 
800 /im, either in SCUBA mapping or UKT14 photometry, or botlfl We have also detected 
17 of these stars at 450/im. Although many of these stars have already been observed 



at millimeter and/or sub- millimeter wavelengths ([Hillenbrand et al.l 1 19921 ; iManningd 11994 



Henning et al.lll998l ; iFuente et al.lll998at iNatta et al 



2004t ?), our study provides the largest 



sub-millimeter sample of HAEBE stars analyzed in a consistent way. We have also included 
Spitzer MIPS photometry for those sources for which such data were publicly available. In 
our analysis we fit a simple isothermal graybody model to the observed SEDs from far- 
infrared to millimeter wavelengths. This allows us to determine the temperature of the cold 
dust, the dust emissivity, and the disk mass. 

We find that disks are ubiquitous in Herbig Ae and F stars, but are very rare in early 
B stars. This does not mean that early B stars do not form with disks, but instead suggests 
that any such disks evolve rapidly and get destroyed quickly by photoevaporation. There are 
currently about five or six early B-stars known to have disks. The three stars we detected in 
this survey — MWC 297, MWC 349, and LkHa 101 — look rather different from the accretion 
disks surrounding Herbig Ae and T Tauri stars. The gas in these disks is largely ionized, and 
the emission at millimeter and sub-millimeter wavelengths is dominated by a thermal wind 
powered by photoevaporation. Only LkHa 101 appears to have a dust excess at millimeter 
wavelengths, suggesting that it is perhaps younger than the other two and has not yet 
had time to photoevaporate the disk. All three stars may have some neutral gas in the 
midplanes of their disks, although there are no detections of molecular gas toward any of 
these three stars at radio wavelengths. MWC 137 and HD 200775, which we did not detect 
in our SCUBA survey, are even fainter and possibly more evolved. Our survey also includes 
two deeply embedded objects (AFGL961A and V1318 Cyg-mm), which have luminosities 
comparable to early B-stars. The disks/envelopes around these stars look very similar to 
Class or Class I protostars, although they appear more massive than disks around young 
low mass stars. Our results therefore suggest that at least some B stars form with accretion 



We note, however, that this is not an unbiased survey, most targets were chosen, because they had a 
high likelihood to be surrounded by accretion disks. 
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disks in an analogous way to how low mass stars form with disks. 

Disks around HAEBE stars appear both larger and more massive than disks around T 
Tauri stars. We are not aware of any accretion disk around a T Tauri star that has been 
resolved with a single dish telescope in sub-millimeter continuum emission. Yet at least 
three HAEBE disks, those around MWC480, HD 163296 and PV Cep, have been resolved 
with SCUBA. The sizes that we find agree well with sizes derived from millimeter-aperture 
synthesis observations. Overall, however, there is no clear correlation between disk size and 
spectral type in our sample, or even any correlation between disk mass and spectral type. 

What we do find, however, is a clear correlation between dust emissivity and disk mass 
in the sense that a lower dust emissivity corresponds to a lower disk mass. This is shown in 
Figure EU where we have taken disk masses and emissivities from the isothermal graybody 
fits given in Table HI The only data we omitted from this plot were the two protostellar 
sources (AFGL961A and V1318 Cyg-mm) and the very distant star V645 Cyg, because 
here the masses are dominated by the envelopes, rather than the accretion disks. As we can 
see from Figure 1541 the correlation is rather good (r = 0.8). What this suggest is that when 
(3 decreases, the mass of the cold dust that we can detect decreases. Note that our modeling 
assumes that the dust opacity is constant at 250 fim, i.e., the "standard" Hildebrand opacity, 
while newer studie s show that lar ge grains, corresponding to low (3- values, have lower opacity 



at 250 //m, see e.g. iDraind (120061 ) . Our disk mass estimates may therefore be underestimated 
for low /3-values. Comparison with other studies, including detailed modeling of some of the 
disks in our sample, indicate that the agreement in disk mass is surpr isingly good, and 



usually within a facto r of two. Since a low f3- value suggests grain growth ( jNatta et al.ll2004 



Rodmann et al.l 120061 ) . the disk mass might not actually be dropping; rather, the dust mass 
could be getting buried into large grains, pebbles, and planetesimals. Once buried, we cannot 
detect it anymore. Thus, low (3, which is correlated with low "detectable" disk mass, might 
be a good indicator of planetesimal growth rather than of decreasing disk mass of solids. All 
stars that show low /3-values are isolated HAEBE stars, i.e., they are relatively old, 3-15 
Myr, which seems to be the time required for significant grain growth. 

A few of these stars have SEDs that suggest they are transition stars, i.e., they show 
weak near-infrared excesses but str ong mid-to-far-infrared excesses, s uggesting that they have 



developed gaps largely free of dust ( iNajita. Strom fc Muzerolldl20071 ). Another characteristic 



feature of transition sources, which has become apparent in this study, is that the millimeter 
and far-infrared regions of the SEDs for these HAEBEs show that the outer portions of the 
disks around these stars have both cold and warm dust components. All transition stars have 
cold outer disks, which dominate the emission in the millimeter and sub-millimeter regime, 
while most of the emission at 60 fim or shorter wavelengths comes from warmer dust (T^ > 
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100 K) in the disk, i .e., presumably emission from the hot surface layers of the disk or the 
inner rim or "wall" (jCalvet et al.l l2002t llsella. Carpenter &: Sargentl 120091 ) . The same trend 
is seen for most of the isolated HAEBE stars, i.e., there is a clear excess at 60 fim requiring 
a warmer dust component in addition to the cold outer disk, which dominates the emission 
at wavelengths longer than 100 ^/m. This is not seen in younger HAEBE disks, like R Mq n 
, PV Cep or AB Aur, or for ex. young FU Ori disks (ISandell. Jessop and Jennessl l200ll ). 
which are all well fit with a single temperature isothermal fit down to 60 //m. A warm dust 
component, producing a clear excess at 60 fim is therefore a clear sign of a disk approaching 
the transition disk stage. 

The strong correlation between disk mass and dust emissivity and lack of correlation 
between spectral type and disk mass suggests that the evolutionary age of the disk affects the 
measurable disk mass much more strongly than does the mass of the central star. HAEBE 
stars do have somewhat larger and more ma ssive disks than T Tauri star s. For example, in 
the large SCUBA survey of Taurus- Auriga by lAndrews fc Williams! ( 120051 ) . the most massive 
T Tauri disk is 0.2 M Q , while there are at least three HAEBE stars in our sample (Table H]) 
that have an equally large or more massive disk. 



8. Summary 

We have studied a large sample of HAEBE stars in the sub-millimeter using SCUBA 
on JCMT. We find that disks around Ae and F stars ubiquitous, while disks around early 
B stars are rare. We have, however, identified several examples of disks in early B stars. 
For very young embedded stars with the luminosity of an early B star, they look similar 
to other young disks, albeit more massive. In early B-stars, which already illuminate H II 
regions, such disks are largely ionized and have very little cold dust. The emission in the 
sub-millimeter in all of these cases is dominated by a thermal wind, presumably driven by 
photoionization of the disk. 

We find several examples of HAEBE stars (LkHa 215, HD 259431, CD-42 1172, MWC 297, 
and TY CrA) for which our imaging survey clearly shows that the emission previously iden- 
tified as dust emission from circumstellar disks around these stars instead originates from 
ionized gas (as in the case of MWC 297), or from nearby younger objects, or from the 
surrounding molecular clouds in which these stars reside. These stars are, in fact, free of 
detectable dust emission at long wavelengths; we conclude that previous SED models of 
these stars that derive properties of their circumstellar disks are incorrect. 

Our survey includes several examples of isolated HAEBEs, which have ages of 3 - 15 
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Myrs and still have gas-rich disks. A few of these are transitional stars, which have already 
formed large inner gaps (perhaps due to planet formation) and show evidence for grain 
growth. Most of the isolated HAEBEs show signs of disk evolution. They have cold outer 
disks, which dominate the emission in the far-infrared and sub-millimeter, but show strong 
excesses at 60 /im due to warm gas in their disks, i.e., they cannot be fitted with a single 
temperature isothermal dust model. 

This work made extensive use of the SIMBAD Astronomical Database at the Centre 
de Donnees astronomiques de Strasbourg, France, and NASA's Astrophysics Data System 
Abstract Service. This research used the facilities of the Canadian Astronomy Data Centre 
operated by the National Research Council of Canada with the support of the Canadian 
Space Agency. This work is based in part on observations made with the Spitzer Space Tele- 
scope, which is operated by the Jet Propulsion Laboratory, California Institute of Technology 
under a contract with NASA. We thank Meredith Hughes for giving us flux densities she 
observed with the SMA and we also acknowledge useful discussions and help from William 
Vacca and Bhaswati Mookerjea. Comments and suggestions by the referee were also appre- 
ciated. 
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Table 1. Stellar parameters of our target stars 
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200 


5,8 


m97bc36 


MWC 137 


V1308 Ori 


06 18 45.46 


+15 16 53.4 


BOep 


3 


1.5xl0 3 


36 


1300 


2 


m99an21,m01bi09 


LkHa 215 


V699 Mon 


06 32 41.80 


+10 09 33.6 


B6-8 


1, 37 


2.95X10 3 


13 


800 


38 


m99bu03 


VYMon 


HBC 202 


06 31 06.94 


+10 26 05.0 


B8:c 


3 


870 


39 


800 


38 


UKT14 data 


HD 259431 


MWC 147 


06 33 05.19 


+10 19 20.0 


B6 


1 


1.07X10 4 


13 


800 


38 


m99bu03 


AFGL 961 A 


IRAS 06319+0415 


06 34 37.6 


+04 12 43.9 


B2-B3/B5 


6 


3xl0 3 


16 


1500 


17 


m98aill 


RMon 


MWC 151 


06 39 09.95 


+08 44 10.8 


B8 Illev 


4 


740 


13 


800 


52,53 


UKT14 data 


LkHa 25 


V590 Mon 


06 40 44.64 


+09 48 02.1 


B7 


1 


63 


13 


800 


13,24 


m00bc32 


HD 135344B 


SAO 206462 


15 15 48.44 


-37 09 15.9 


F4Ve 


7 


8.5 


13 


140 


7, 47 


m96bu90 


HD 141569 


SAO 140789 


15 49 57.75 


-03 55 16.4 


AO Vc 


4 


23 


13 


99 


13,5 


m00ua06 


HD 142666 


V1026 Sco 


15 56 40.02 


-22 01 40.0 


A8 Vc 


25 


17 


13 


145 


13,25 


m97bc36 


HD 144432 


SAO 184124 


16 06 57.96 


-27 43 09.8 


A9 IIIc 


8 


15 


13 


~145 


13,25 


m97bc36 


HR5999 


HD 144668 


16 08 34.29 


-39 06 18.3 


A6 Ille 


8 


115 


13 


208 


13 


m01bi09 


HD 150193 


MWC 863 


16 40 17.92 


-23 53 45.2 


A3 IIIc 


8 


49 


13 


150 


13 


m01bi09 


CD-42° 11721 


V921 Sco 


16 59 06.77 


-42 42 08.4 


BOep 


40,41 


2.34X10 4 


16 


1150 


42 


m99au36 


KK Oph 


HBC 273 


17 10 08.07 


-27 15 18.2 


A5 Ve 


8 


26 


13 


160 


13,26 


m01bi09 


HD 316285 


MWC 272 


17 48 14.04 


-28 00 53.2 


Be 


43 


2.8 xlO 5 


43 


1850 


43 


m98au64 


HD 163296 


MWC 275 


17 56 21.29 


-21 57 21.9 


A3Vc 


25,50 


35 


49 


122 


5 


m97bn2 1 ,m0Xa /bec05 


HD 169142 


MWC 925 


18 24 29.78 


-29 46 49.4 


A7 Vc 


8 


22 


8,49 


145 


25,47 


m97bc36 


MWC 297 


NZ Ser 


18 27 39.53 


-03 49 52.2 


B1.5 Vc 


4 


1.07X10 4 


19 


250 


59 


m97bi09 


MWC 300 


V431 Set 


18 29 25.69 


-06 04 37.3 


B[e] 


3 


1.3xl0 5 


20 


1800 


20 


m99au36 


HD 176386 


CD-37° 13023 


19 01 38.93 


-36 53 26.6 


A0V 


48 


36 


48 


130 


47 


m00au58 


TY Cr A 


CD-37° 13024 


19 01 40.83 


-36 52 33.9 


B8V 


45,46 


67 


46 


130 


47 


m00au58 


BD+40°4124 


V1685 Cyg 


20 20 28.37 


+41 21 51.4 


B3 


1 


6.0X10 4 


13 


980 


5 


m97bn21 


V 1686 Cyg 


LkHa 224 


20 20 29.26 


+41 21 28.6 


F9 


1 


257 


13 


980 


5 


rri97bn21 


V1318CygS 


LkHa 225 S 


20 20 30.59 


+41 21 26.0 


A5-Fe 


2 


1600 


22 


980 


5 


m97bn21 



Table 1 — Continued 



Star 


Alt. 


a(2000.0) 


5(2000.0) 


Spectral 


rcf. 




ref. 


Distance 


rcf. 


Project ID 




Name 


[h m s] 


[° ' "1 


Type 




[L©] 




[pc] 






MWC 349 


V 1478 Cyg 


20 32 45.53 


+40 39 36.6 


Be 


44 


3xl0 4 


44 


1200 


44 


m01bi09 


PV Cep 


HBC 696 


20 45 53.96 


+67 57 38.9 


A5 


1 


80 


54 


500 


56 


m97bn21, & UKT14 data 


HD 200775 


MWC 361 


21 01 36.92 


+68 09 47.8 


B3 


1 


1.9 x 10 3 


36 


430 


5 


m99bu03,m01bi09 


V 645 Cyg 


RAFGL 2789 


21 39 58.24 


+50 14 21.2 


Bl-2 


57 


2 - 6xl0 4 


57 


4200 


57 


m01bi09 


MWC 1080 A 


V628 Cas 


23 17 25.59 


+60 50 43.6 


BOeq 


30 


4.6xl0 4 


54 


2200 


55 


m01bi09 



(1) Hernandez et al. (2004), (2) Hillenbrand et al.(1992), (3) The et al. (1994), (4) Mora et al. (2001), (5) van den Ancker et al. (1998), (6) Castelaz ct al. (1985), (7) van Boekel et al. 
(1995), (8) Blondcl et al. (2006), (9) Chavarria-K (1985), (10) Natta et al. (1992), (11) Monnier et al. (2008), (12) Torres et al. (2009), (13) Manoj et al. (2006), (14) Tuthill et al. (2002), 
(15) Bertout et al. (1999), (16) this paper (17) Williams et al. (2009), (18) Habart et al. (2003), (19) Alonso-Albi et al. (2009), (20) Miroshnichenko ct al. (2004), (21) Siebenmorgen et al. 
(2000), (22) Marvel (2005), (23) Reipurth et al. (1997), (24) Finkenzellcr & Mundt (1984), (25) van Boekel et al. (2005), (26) Leincrt ct al. (2004), (27) Lagage et al. (1993), (28) Natta et 
al. (1992), (29) Strom & Strom (1994), (30) Cohen & Kuhi (1979), (31) Herbig ct al. (2004), (32) Mccus ct al. (2001), (33) Dent et al. (2005), (34) Merin ct al. (2004), (35) Acke et al. 

(2004), (36) Bcrilli ct al. (1992), (37) Valcnti ct al. (2003), (38) Hcrbst ct al. (1982), (39) Casey & Harper (1990), (40) do Winter & The (1990), (41) Shore ct al. (1990), (42) Borgcs 
Fernandcs et al. (2007), (43) Hillicr ct al. (1998), (44) Cohen ct al. (1985), (45) Casey et al. (1995), (46) Casey ct al. (1998), (47) dc Zceuw ct al. (1999), (48) Meyer & Wilking (2009), 
(49) Ackc et al. (2004), (50) Gray & Corbally (1998), (51) Cohen, Harvey & Schwartz (1985), (52) Walker (1956), (53) Close et al. (1997), (54) Evan et al. (1986) (55) Canto et al. (1994), 

(56) Cohen et al. (1981), (57) Miroshnichenko ct al. (2009), (58) Duncin & Crawford (1998), (59) Drew et al. (1997) 
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Table 2a. Positions, deconvolved sizes and integrated flux densities of target stars 



HAEBE 


a(2000.0) 


5(2000.0) 


On X Oh 


P.A. 


S(850 (im) 


S(450 (im) 


star 


[h Ml s] 


ro / Mi 
L J 




L J 


[m.Tvl 
[iiiu y j 


[m.Tvl 
[iiiu y j 


LkHa 198 a 


00 11 26.23 


+58 49 32.2 


11'.' 4 x 5'.' 6 


+29 ± 10 


220 ± 30 




V 376 Cas b 


00 11 26.36 


+58 50 01.4 


< 6" x 6" 




140 ± 20 




V 892 Tau 


04 18 40.64 


+28 19 16.1 


point source 




630 ± 20 


2250 ± 100 


AB Aur 


04 55 45.89 


+30 33 05.1 


point source 




350 ± 20 


2200 ± 100 


MWC 480 


04 58 46.45 


+29 50 39.0 


6'.' 4 x 4'.' 7 


+138 ± 4 


780 ± 20 


3300 ± 290 


HD 34282 


05 16 00.48 


-9 48 35.3 


point source 




360 ± 9 


2700 ± 400 


HD 36112 


05 30 27.56 


+25 19 56.9 


point source 




197 ± 10 


1540 ± 50 


HD 35187 


05 24 01.17 


+24 57 37.6 


point source 




70.6 ± 6 


340 ± 85 


LkHa 101 c 


04 30 14.38 


+35 16 23.6 


6'.' 5 X 2'.' 6 


+26 ± 20 


818 ± 50 




AFGL 961 A 


06 34 37.54 


+04 12 44.6 


12'.' 2 X 7'.'0 


-41 ± 2 


1950 ± 60 


14900 ± 300 


AFGL961 SMM 


06 34 36.71 


+04 12 50.9 


8'.'8 x 5'.'6 


-5 + 10 


1330 ± 50 


12300 ± 300 


HD 135344 B 


15 15 48.56 


-37 09 15.8 


point source 




490 ± 10 


3180 ± 230 


HD 141569 


15 49 57.75 


-03 55 16.4 






10.9 ± 1.3 


64.9 ± 13.3 


HD 142666 


15 56 40.01 


-22 01 39.0 


< 2'.' 2 




313 ± 5 


1140 ± 35 


HD144432 


16 06 57.90 


-27 43 09.8 






129 ± 6 




HR5999 


16 08 34.42 


-39 06 18.3 


point source 




64 ± 20 




HD 150193 


16 40 17.85 


-23 53 44.8 


point source 




101 ± 10 




KK Oph 


17 10 08.12 


-27 15 18.8 


point source 




91 ± 15 




HD 163296 


17 56 21.42 


-21 57 21.6 


4'.' 6 X 3'.' 2 


-13 ± 27 


1900 ± 200 


8700 ± 370 


HD 169142 


18 24 29.87 


-29 46 48.9 


~ 2'.'3 x 1'.'6 




565 ± 10 


3340 ± 115 


MWC 297 


18 27 39.53 


-03 49 50.4 


<4" 




570 ± 30 


1450 ± 200 


V1318 CygS-mm 


20 20 30.67 


+41 21 29.9 


8" x 3" 


+62 ± 1 


2050 ± 50 


10200 ± 200 


MWC 349 


20 32 45.53 


+40 39 36.6 


point source 




2600 ± 70 


5000 ± 1100 


PVCep 


20 45 54.31 


+67 57 40.7 


4'.'9 x 2'.'8 


+74 ± 10 


1000 ± 20 


6550 ± 150 


V645 Cyg 


21 39 58.26 


+50 14 22.3 


14'.' 2 x 8'.' 5 


-65 ± 3 


2400 ± 60 


20800 ± 600 


MWC 1080 A 


23 17 24.35 


+60 50 40.4 


point source 




250 ± 50 





a LkHa 198 and LkHa 198 B are blended with the nearby Class source, LkHa 198-mm, resulting in uncertain 
size and position. LkHa 198 B appears to dominate the sub-millimeter emission. 

b Peak offset by 3" from the optical position of V376Cas, possibly emission from the surrounding cloud. 

c Flux density at 750 fan, S(750^m) = 940 ± 100 mjy 
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Table 2b. Submillimeter photometry of compact HAEBE stars 



HAEBE 


A 


HPBW 


Flux density 


Date 


star 


H 


["] 


[mjy] 


mm/dd/yr 


AB Aur 


1300 a 


32.0 


136 ± 15 


12/10/89 




1100 


18.4 


192 ± 19 


08/20/88 




800 


15.8 


534 ± 69 


08/20/88 




800 


16.8 


475 ± 56 


08/16/90 




450 


17.5 


1760 ± 360 


08/16/90 


R Mon 


1300 a 


32.0 


96 ± 15 


12/09/89 




1300 


19.0 


36 ±18 


10/09/91 




1100 


18.4 


80 ± 20 


08/19/88 




850 


17.8 


171 ± 24 


08/16/90 




800 


15.8 


260 ± 30 


08/19/88 




800 


16.8 


197 ± 30 


08/16/90 




450 


17.5 


1200 ± 160 


08/19/88 




350 


19.0 


2600 ± 500 


08/19/90 


VY Mon 


1300 


19.0 


189 ± 23 


10/09/91 




1100 


18.5 


247 ± 27 


10/09/91 




800 


16.5 


552 ± 55 


10/09/91 


KK Oph 


1100 


18.4 


36 ± 11 


08/11/91 




800 


16.8 


109 ± 20 


08/18/90 




800 


16.4 


119 ± 53 


10/06/91 


PV Cep 


1300 


32.0 


0.375 ± 0.044 


12/10/89 




1300 


19.5 


0.36 ± 0.023 


08/18/90 




1100 


18.4 


0.50 ± 0.04 


08/20/88 




1100 


18.4 


0.47 ± 0.03 


10/27/88 




850 


17.8 


1.01 ± 0.07 


08/18/90 




800 


15.8 


1.3 ± 0.1 


08/20/88 




800 


15.8 


1.11 ± 0.09 


10/27/88 




800 


16.8 


1.16 ± 0.06 


08/18/90 




800 


13.5 


1.2 ± 0.1 


08/20/88 




750 


17.5 


1.41 ± 0.09 


08/11/91 




450 


17.5 


7.91 ± 0.62 


08/18/90 




450 


17.5 


3.48 ± 1.20 


08/11/91 




450 


13.5 


10.6 ± 0.8 


08/20/88 




450 


7.9 


10.5 ± 0.5 


08/20/88 




350 


17.7 


17.6 ± 1.0 


08/20/88 




350 


10 


16.2 ± 0.7 


08/20/88 



a CSO observations 
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Table 2c. Non-detections 



Star 


850 [im upper limit 
[mJy/beam] 


Comments 


MWC 137 


15 


cloud emission, nearby sources 


LkHa215 


22 


emission from the surrounding nebula 


HD 259431 


20 


emission from the surrounding nebula 


LkHa 25 


25 


cloud emission 


CoD-42° 11721 


50 


nearby sources 


HD 316285 


70 


cloud emission 


MWC 300 


7 


empty field 


HD 176386 


50 


cloud emission 


TY CrA 


50 


cloud emission 


BD+40°4124 


35 


cloud emission 


V1686 Cyg 


35 


cloud emission, nearby source 


HD 200775 


20 


emission from the surrounding nebula 
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Table 3. Spitzer MIPS flux densities of program stars 



Star 


24 


70 /zm 


AOR ID 




[Jy] 


[Jy] 




AB Aur 


s a 


s a 


rl2662784 


V892 Tau 


s a 


25.41 ± 0.10 b rllll9696, rll229952 


AFGL961 


s a 


S a 


rl2234249, rl2234752 


HR5999 


2.89 ± 0.72 


3.78 ± 0.36 


c2d 


HD 141569 


1.47 ± 0.01 


4.70 ± 0.02 


rlll81312 


HD 150193 


S a 


3.05 ± 0.29 


c2d 


CD-42 11721 


S a 


S a 


MIPSGAL 


KK Oph 


S a 


s a 


rl4094848 


BD+40°4124 


S a 




r22512128, r22512384 


V1686 Cyg 


S a 




r22512128, r22512384 


V1318CygS-mm 


S a 


85.68 ± 0.66 


r22512128, r22512384 


MWC 349 


S a 


6.72 ± 0.30 


r22510336 



a Saturated 
b Mildly saturated 
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Table 4. Isothermal graybody fits to HAEBE "disks" 



Source 


A - a 

[Kl 






Oi 








■ LyL tot 
L 1VJ -0J 


^aust 

[Lnl 
L ©J 


V892 Tau 


40 


2. 


.19 


± 





.09 





.52 


0.009 


0.27 


AB Aur 


66 


2. 


.89 


± 


0. 


.04 


1 


.05 


0.006 


5.0 


MWC480 


28 


2. 


,65 


± 





,06 





,81 


0.02 


0.23 


HD 34282 


29 


3. 


.03 


± 


0. 


.12 


1 


.29 


0.18 


1.85 


HD 35187 


44 


2 


,64 


± 





,34 





,87 


0.002 


0.15 


HD 36112 


43 


3 


.03 


± 


0. 


.10 


1 


.25 


0.014 


1.0 


AFGL961 A 


37 


2 


,73 


± 





,02 


1 


.51 


12.5 


930 


R Mon 


62 


3 


,08 


± 


0. 


,15 


1 


,3 


0.12 


110 


VY Mon a 


38 


2. 


.31 


± 


0. 


.93 


1 


.9 


1.9 


102 


HD 135344 B 


32 


2. 


.94 


± 





.30 


1 


.38 


0.025 


0.20 


HD 141569 


40 


2. 


,21 


± 


0. 


.06 


0. 


.5 


0.00009 


0.03 


HD 142666 


31 


2 


.42 


± 





.07 





.72 


0.005 


0.08 


HD 144432 


40 


2 


,70 


± 


0. 


,39 


0. 


.59 


0.002 


0.10 


HD 150193 


44 


3 


.03 


± 


0. 


.15 





.65 


0.002 


0.12 


KK Oph 


40 


1 


.70 


± 


2. 


.12 


0. 


.74 


0.002 


0.13 


HD 163296 


28 


2 


,66 


± 





,05 





,94 


0.065 


0.52 


HD 169142 


41 


2 


,84 


± 


0. 


.06 


1. 


.30 


0.03 


0.50 


HR5999 


41 


3 


,03 


± 


5 


,18 


1 


.75 


0.009 


1.1 


V645 Cyg 


38 


2 


,84 


± 


0. 


.06 


1. 


.55 


124 


11750 


V1318 Cyg-mm 


41 


2 


,55 


± 


0. 


.06 


0. 


.85 


1.9 


118 


PV Cep 


29 


3. 


.09 


± 





,06 


1 


,59 


1.0 


13.7 



a Fit to UKT14 photometry data, likely to include some con- 
tribution from the surrounding cloud. 
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Fig. 1. — Isothermal graybody fits to SCUBA and literature data for isolated HAEBE stars. 
Flux densities used in the fits are marked by filled triangles, while data points marked with 
open circles were excluded from the fits. In the top left corner of each plot we give the name 
of the star and essential fit results (T d , (5). 
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10 100 1000 10 4 10 100 1000 



v [GHz] v [GHz] 

Fig. 2. — Isothermal graybody fits to isolated HAEBE stars continued. Labeling as in 
previous figure. In the right panel for HD 169142, we show a isothermal fit with two dust 
components at different temperatures. The cold dust, which dominates at the long wave- 
lengths is shown as a dashed line, the warmer dust as a dash-dotted line, and the combined 
two-temperature fit as a solid line. 
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Fig. 3.— Deconvolved sub-millimeter images of MWC480 at 850 and 450 fim plotted 
in gray-scale and overlaid with contours. The contour levels are linear with six contours 
between the lowest contour level (3-cr) and the peak flux density. The 3-a level is at 54 
mJy beam -1 and at 840 mJy beam -1 for 850 and 450 //m, respectively. The disk is resolved 
by SCUBA both at 850 and 450 /im. The position of MWC480 is shown with a star symbol. 
The HPBW is plotted in the bottom left corner of each image. 
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Fig. 4. — Deonvolved sub-millimeter images HD 163296 at 850 and 450 /im. The contours 
are plotted with linear steps. At 850 /im we plot six contours between the lowest contour 
level, 32 mJy beam -1 (2-cr) and the peak flux density, 1.8 Jy beam -1 . At 450 /im we plot six 
contours starting from 0.45 Jy beam -1 to 7.24 Jy beam -1 . The disk is resolved by SCUBA 
both at 850 and 450 /im. The position of HD 163296 is marked by a star symbol. The 
HPBW is plotted in the bottom left corner of each image. 
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Right Ascension (2000.0) 

Fig. 5. — Deonvolved SCUBA images of HD 169142 in grayscale overlaid with contours. The 
contour levels are linear for both 850 and 450 /im with six contours starting at 2-cr. At 850 
yum contours go from 15 mJy beam -1 to 537 mJy beam -1 , at 450 //m from 0.28 Jy beam -1 
to 3.17 Jy beam -1 . The disk is marginally resolved (see text). 
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MWC349 



450 yu,m 




MIPS 70 yum 



0.5 



1.5 



2.5 

log(z/) [GHz] 



3.5 



4.5 



Fig. 6. — Least squares fit to radio and millimeter data for MWC 349 showing that the 
star has a thermal wind spectrum which dominates the emission to about 450 /im with the 
flux density proportional to the frequency z/°- 67±0 03 . If we assume that the wind is launched 
at 15 AU from the star, the wind should become optically thin at 450 ^m (see text). 
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Fig. 7. — SCUBA 850 //m image in grayscale overlaid with contours of the dense cloud in 
which LkHa 198 and V376 Cas are embedded. The cloud core is quite extended and probably 
more extended to the northeast, where the image goes negative because we chopped on 
extended emission. The emission is dominated by a Class source, LkHa 198-mm, which does 
not have an optical or near-infrared counterpart. This image suggests that both LkHa 198, 
its infrared companion and V376 Cas may be associated with dust emission, but we do not 
have high enough spatial resolution to separate the stars from the surrounding strong cloud 
emission. The contours are plotted with six logarithmic intervals from 40 mJy beam" 1 to 
750 mJy beam -1 . The HPBW is shown in the bottom right corner. 
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04 h 18 m 44 s 43 s 42 s 41 s 40 s 39 s 38 s 37 
Right Ascension (2000.0) 



Fig. 8. — Deconvolved SCUBA images of V 892 Tau in grayscale overlaid with contours. The 
contour levels are evenly spaced with six contours between the lowest contour level, ~ 1-a 
and the peak flux density. At 850 /xm contours go from 36 mJy beam -1 to 630 mJy beam -1 , 
at 450 /iin from 0.2 Jy beam -1 to 2.25 Jy beam -1 . The position of V892 Tau is marked by 
a star symbol. The HPBW is shown in the bottom right corner. 
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04 h 18 m 55 s 50 s 45 s 40 s 35 s 30 s 



Right Ascension (2000.0) 

Fig. 9. — MIPS 70 /iin image of V892 Tau in grayscale overlaid with logarithmic contours 
and covering a much larger field than the SCUBA images. The extent of the SCUBA images 
is drawn as a dotted rectangle. At 70 /iin there is extended nebulosity east and northwest 
of the star, but no sign of embedded sources in the nebulosity. Two faint 70 fim. sources are 
visible within the displayed area. The 70 ^m HPBW is shown in the bottom left corner. 
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Fig. 10. — Isothermal graybody fit of V892 Tau assuming a dust temperature of 40 K. 
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Fig. 11. — SCUBA 850 /im grayscale image of LkHa 101. The contours are linear with the 
lowest contour at 100 mJy beam -1 and a step of 100 mJy beam -1 . The peak flux is 0.9 
Jy beam -1 . The position of LkHa 101 is marked on the figure. The HPBW is shown in the 
bottom right corner. 
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Fig. 12. — Power-law fit to VLA A-array data for LkHa 101 predicting a spectral index, a 
= 0.71 ± 0.07 for the thermal wind. If we assume that the wind is launched at 10 AU from 
the star, the wind will become optically thin at ~ 150 GHz. 
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Fig. 13. — Deconvolved sub-millimeter images of AB Aur at 850 and 450 ^m plotted in 
gray-scale and overlaid with contours. The disk is unresolved both at 850 and 450 ^m. The 
peak flux density is 0.349 and 2.2 Jy beam" 1 for HPBWs of 14" and 8" for the 850 and 450 
fim images, respectively. The contour levels are logarithmic with six contours between the 
2a noise level and peak flux density. The 2a noise level is 18 mJy beam -1 at 850 /im and 
110 mJy beam -1 at 450 /J,m. The position of AB Aur is marked by a star symbol. The faint 
extension to the southwest in the 450 image appears to be spurious, because it does not 
agee in position with the faint extension seen at 850 /im. 
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Fig. 14.- 



An isothermal graybody fit to the SED for AB Aur. 
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Fig. 15. — CO J = 2 — > 1 spectrum of VY Mon showing high velocity emission from about 
-28 to +18 km s _1 . 
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Fig. 16. — Map of CO J = 2 — > 1 high velocity emission around VY Mon. The solid contours 
show the red-shifted emission integrated over the velocity range —14 to —3 km s -1 with the 
lowest contour at 2 K km s" 1 and steps of 6 K km s -1 . The dashed contours show the 
blue-shifted emission integrated over the velocity range +1 to +12 km s -1 with the lowest 
contour at 3 K km s _1 and steps of 5 K km s _1 . The bipolar outflow appears unresolved 
and centered on VY Mon, which is marked by a star symbol. 
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Fig. 17. — Upper panels: Contour plots of 850 and 450 fim SCUBA images. The 850 /zm 
image (left) is overlaid on a MIPS 70 jum image in grayscale, also shown below. The triangles 
in the saturated region are AFGL961A (furthest to the left) and AFGL961B (second from 
left). The other two triang l es are SMA3 (second from right) and source C (furthest to the 
right) from IWilliams et al.l ( 120091 ) . At 850 /im we plot ten evenly spaced contours from 0.15 
Jy beam -1 to 2.3 Jy beam -1 . At 450 /im we plot eight evenly spaced contours between 1. 
Jy beam -1 and 10 Jy beam -1 . Lower panels: MIPS 70 and 24 /im images in grayscale. 
Both images are severely saturated (pure white in grayscale). At 24 /im source AFGL961 C 
breaks up into a double source. The source immediately to the southeast of C is saturated 
in the 24 /im MIPS image. The HPBWs are shown in the bottom left corner of each image. 
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Isothermal graybody fit of AFGL 961 A. 



-81 - 



10 



15 



10 14 r 



10 



10 



13 



12 



in 



Ho 10 



10 a - 



10 c 



10' 



10 c 



R Mon T d = 62 K, |8 = 1.30 



CD 



6 [ g , , , , 




10 



_l I I I I I I L 



J I I I I I I L_ 



100 r n 1000 

v [GHz] 



10" 



Fig. 19. — Isothermal graybody fit of R Mon to free-free correct ed (sub-) millimeter data and 
KAO far infrared data. The 1.3 mm flux density observed by (iFuente et all 120061 ) appears 
anomalously low and was not used in the fit. 
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Fig. 20.— SCUBA 850 and 450 fxm images in grayscale overlaid with contours of the 
MWC 297 region. At 850 /xm we plot ten evenly spaced contours between 30 mJy beam -1 
and 611 mJy beam" 1 , at 450 fim six evenly spaced contours between 0.2 Jy beam" 1 and 
2.0 Jy beam" 1 . MWC 297 is marked by a star symbol. The HPBW is shown in the bottom 
right corner of each image. 



- 83 - 



CO 

cn 

O 



MWC297 



450/^m 



850^171 



115 GHz 




1.5 2 
logO) [GHz] 

Fig. 21.— Least-squares fit to VLA and BIMA data (open squares) up to 3 mm (see 
text), predicting a spectral in dex a ~ 1. In this fit we omitted the VLA data point at 43 
GHz (lAlonso-Albi et al.ll2009l ). which appears anomalously low, suggesting that a significant 
fraction of the emission at this wavelength has been resolved out by the interferometer. 
There may be some dust excess at 850 /im and almost certainly at 450 jttm. The dashed line 
shows that the free-free emission will flatten out when the whole jet becomes optically thin. 
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Fig. 22. — SCUBA 850 and 450 /xm images of the BD 40°4124 field in grayscale overlaid with 
contours. The contour levels are evenly spaced with 10 contours from the 3-a level to the 
peak flux density. At 850 fim the lowest contour is at 40 mJy beam" 1 and the highest at 2.45 
Jy beam -1 . At 450 fxm the lowest contour is at 300 mJy beam -1 with the highest contour at 
8.9 Jy beam -1 . The stars, BD 40°4124, V1686 Cyg and V 1318 Cyg N and S are marked by 
star symbols and the first two are labeled on the 450 /iin image. V 1318 Cyg N and S coincide 
with a strong sub-millimeter source, although the emission is likely completely dominated 
by an invisible Class source situated about one arcsecond northeast of V 1318 Cyg S (see 
text). The HPBWs are shown in the bottom right corner of each image. 
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Fig. 23. — SCUBA 850 and 450 fj,m images of PV Cep in grayscale overlaid with contours. 
The contours are logarithmic with the lowest contour at the 3-a level and the highest at the 
peak flux density. At 850 fim the lowest contour is at 50 mJy beam -1 with the peak flux 
density at 1 Jy beam -1 . At 450 /iin the corresponding values for the contour levels are 300 
mJy beam -1 and 5.5 Jy beam -1 . PV Cep is marked by a star symbol and H 2 masers with 
triangles. The circumstellar disk/envelope is well resolved in these SCUBA images with a 
size of ~3'.'6. We also see faint emission from the surrounding dust cloud in which PV Cep 
is embedded. The HPBWs are shown in the bottom right corner of each image. 



-86 - 




24.- 



An isothermal graybody fit to the SED for PV Cep. 
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Fig. 25. — SCUBA 850 /im image of V645 Cyg in grayscale overlaid with contours. The 
contours are logarithmic with eight contours between 40 mJy beam -1 and 1.78 Jy beam -1 . 
The HPBW is shown in the bottom right corner. 
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Fig. 26. — SCUBA 850 /im image of the MWC 1080 field in grayscale overlaid with con- 
tours. The contours are logarithmic with si x contours between 150 mJy beam -1 and 700 
mJy beam -1 . The VLA sources 3, 4 and 5 ( IGirart et al.l 120021 ) are marked with triangles. 
VLA 4 coincides with the B0 star MWC 1080 A, also marked with a star symbol. The HPBW 
is shown in the bottom right corner. 
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Fig. 27. — An 850 /zm contour map of MWC 137 overlaid on an optical image in grayscale. 
The contour levels are linear starting at 50 mJy beam -1 with a step of 50 mJy beam -1 . 
MWC 137 is marked by a star symbol. The HPBW is shown in the bottom right corner. 
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Fig. 28. — An 850 /zm contour map of the reflection nebula NGC 2245 overlaid on an optical 
image in grayscale. The contour levels are linear starting at 60 mJy beam -1 with a step of 
60 mJy beam -1 . LkHa215 is not detected at 850 fim. Northwest of the star and outside 
the reflection nebula there is a compact sub-millimeter source with an 850 flux density 
of 150 mJy, which has neither an optical nor a near-infrared counterpart. The HPBW is 
shown in the bottom right corner. 
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Fig. 29. — An 850 /zm contour map of the reflection nebula NGC 2247 overlaid on an optical 
image in grayscale. The contour levels are linear, starting at 50 mJy beam -1 with a step of 
50 mJy beam -1 . The nebula is illuminated by HD 259431, which is close to a minimum in 
the dust emission. Two faint sub-millimeter sources, marked with black plus signs, are seen 
outside the nebula. The HPBW is shown in the bottom right corner. 
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Fig. 30.— Upper panels: 850 and 450 SCUBA maps of the CD-42° 1172 field. The 
position of CD— 42° 1172 is marked by a star symbol, and the two 13 cm radio sources 
( jThompson. Urguhart fc White! 120041 ) are marked by filled squares. The contour levels are 
logarithmic. At 850 /iin there are seven contour levels between 0.1 Jy beam" 1 and 1.58 
Jy beam -1 , at 450 fx they go from 1 Jy beam -1 to 15.8 Jy beam -1 The SCUBA HPBWs 
are shown in the bottom right corner of each panel. Lower panels: False color images of 
24 and 70 /im MIPS images overlaid with the same 850 /im contours as shown in the left 
upper panel. Note that the 24 /im image is saturated on CD— 42° 1172 and the nebulosity 
surrounding it. The source to the northwest of the star is also saturated. 
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Fig. 31. — Sub-image of large SCUBA 850 /im scan map showing HD 176386 and TY CrA. 
Although both stars are definitely embedded in a dense dust cloud, neither of them are 
associated with any sub-millimeter continuum emission. The HPBW is shown in the bottom 
right corner. 
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Fig. 32. — Left panel: Large 850 /mi scan map of the reflection nebula NGC 7023 plotted 
with contours and overlaid on an optical image in gray scale. We applied a logarithmic 
stretch to the optical image to enhance the faint nebulosity. The contour levels are linear, 
starting at 40 mJy beam -1 with ten contours up to the peak flux density of 712 mJy beam -1 . 
The HPBW is shown in the bottom right corner. Right panel: A much deeper 850 //m 
jiggle map with a 3-cr rms of 25 mJy beam -1 . The lowest contour in this map is at 25 
mJy beam -1 . The rest of the contours are plotted logarithmically to the peak flux density, 
530 mJy beam -1 . Even though HD 200775 falls inside a 3-cr contour, we consider it a non- 
detection, because it is not distingui shable from other nois e features in the map. However, 
our observations are consistent with lOkamoto et al.l (120091 ). who detected HD 200775 with 
the SMA at 350 GHz (856 fim) with a flux density of 35 ± 5 mJy, The SCUBA HPBW is 
shown in the bottom right corner of the image. 
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Fig. 33. — Large SCUBA scan maps at 850 and 450 in the direction of HD 316285. 
There is a faint peak at the position of the star at 850 /im (hidden by the star symbol), but 
not at 450 /im. The contours are logarithmic for both images. For the 850 /im image, we 
plot seven contour levels between the lowest contour, 200 mJy beam -1 , and the peak flux 
density, 3.16 Jy beam -1 . For 450 /im we plot eight contours between 1 Jy beam -1 and 25.1 
Jy beam -1 . The HPBWs are shown in the bottom right corner of each image. 
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Fig. 34. — Disk mass as a function of dust emissivity, j3. The error bars show a worst case 
uncertainty of a factor of two for disk mass. The errors for (3 are of the order of 0.1 - 0.2. 
The two extreme outliers, HD 141569 and HR5999 are labeled in the Figure. 



